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fc-EXTENDIBILITY OF HIGH-DIMENSIONAL BIPARTITE QUANTUM STATES 


CECILIA LANCIEN 


Abstract. The idea of detecting the entanglement of a given bipartite state by searching for symmetric extensions 
of this state was first proposed by Doherty, Parrilo and Spedialeri. The complete family of separability tests it 
generates, often referred to as the hierarchy of k-extendibility tests, has already proved to be most promising. The 
goal of this paper is to try and quantify the efficiency of this separability criterion in typical scenarios. For that, 
we essentially take two approaches. First, we compute the average width of the set of /c-extendible states, in order 
to see how it scales with the one of separable states. And second, we characterize when random-induced states are, 
depending on the ancilla dimension, with high probability violating or not the fc-extendibility test, and compare 
the obtained result with the corresponding one for entanglement vs separability. The main results can be precisely 
phrased as follows: on 0 C^, when d grows, the average width of the set of /c-extendible states is equivalent 
to {2/^/k)/d, while random states obtained as partial traces over an environment of uniformly distributed pure 
states are violating the fc-extendibility test with probability going to 1 if s < ((fc — l)^/4/c)d^. Both statements 
converge to the conclusion that, if k is fixed, /c-extendibility is asymptotically a weak approximation of separability, 
even though any of the other well-studied separability relaxations is outperformed by A:-extendibility as soon as k is 
above a certain (dimension independent) value. 


1. Introduction 


Deciding whether a given bipartite quantum state is entangled or separable (or even just close to separable) 
is known to be a computationally hard task (see m and HZ]). Several much more easily checkable necessary 
conditions for separability do exist though, the most famous and widely used ones being perhaps the positivity of 
partial transpose criterion m, the realignment criterion m or the A:-extendibility criterion |15| . All of them have 
in common that verifying if a given state fulfils them or not may be cast as a Semi-Definite Program (SDP) and 
hence be efficiently solved (see e.g. the quite extensive review [U] for much more on that topic). 

We focus here on a relaxation of the notion of separability of quite different kind: the so-called fc-extendibility 
criterion for separability, which was introduced in m- It is especially appealing because it provides a hierarchy 
of increasingly powerful separability tests (expressible as SDPs of increasing dimension), which is additionally 
complete, meaning that any entangled state is guaranteed to fail a test after some finite number of steps in the 
hierarchy. Let us be more precise. 


Definition 1.1. Let fc G N. A state pab on a bipartite Hilbert space A 0 B js k-extendible with respect to B if 
there exists a state Pab'" on A 0 B®*^ which is invariant under any permutation of the B subsystems and such that 
PAB = TrBfe-i Pab*" • 

Theorem 1.2 (The complete family of /c-extendibility criteria for separability, [H]). A state on a bipartite Hilbert 
space A 0 B is separable if and only if it is k-extendible with respect to B for all /c S N. 


Note that one direction in Theorem |1.2| is obvious, namely that a separable state on some bipartite system is 
necessarily /c-extendible for all /c S N (with respect to both subsystems). Indeed, if p = ® Er is separable, 

then ^Tx and J^xPxO’x^'rf’^ are symmetric extensions of p to fc copies of the first and second subsystems 


respectively. The other direction in Theorem 1.2 follows from the quantum finite de Finetti theorem (see e.g. |24I[T^ 
for the seminal statements). The latter establishes, roughly speaking, that starting from a permutation-invariant 
state on some tensor power system and tracing out all except a few of the subsystems, one gets a state that may 
be well-approximated by a convex combination of tensor power states (with a vanishing error as the initial number 
of subsystems increases). 

It is easy to see that if a state is /c-extendible for some /c € N, then it is automatically /c'-extendible for all k' ^ k. 
Hence, the necessary and sufficient condition for separability provided by Theorem |1.2| actually decomposes into a 
series of increasingly constraining necessary conditions for separability, which are only asymptotically also sufficient 
(see Figure]^. In real life however, checks can only be done up to a finite level in this hierarchy. It thus makes 
sense to ask, given a finite /c S N, how “powerful” the /c-extendibility test is to detect entanglement. 

Actually, various more quantitative versions of Theorem |1.2| do exist, that put bounds on how far a /c-extendible 
state can be from separable. Let us mention two quite different statements in that direction. The original result, 
appearing in [T^, establishes that a state on A 0 B which is /c-extendible with respect to B is at distance at most 
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Figure 1. The nested and converging sequence of fc-extendibility relaxations of separability, fc € N 



2d‘^/k, in 1-norm, from the set of separable states. It is a direct consequence of one of the quantitative versions 
of the quantum finite De Finetti theorem. A more recent result, proved essentially in and improved in m, 
stipulates that such a state is at distance at most -^2 in LOCC^ -norm, from the set of separable states 

(see e.g. [3] for a precise definition of the operational one-way-LOCC norm). It relies on the observation that a 
fc-extendible state has a small squashed entanglement, and therefore cannot be distinguished well from a separable 
state by local observers. The main problem of such estimates is that they become non-trivial only when fc 3> de or 
A: ^ In dA- So in the case where dx, de are “big”, can anything interesting still be said for a “not too big” fc? On the 
other hand, these bounds valid for any fc-extendible state are known to be close from optimal (there are examples 
of fc-extendible states whose closest separable state is at distance of order de/fc in 1-norm or of order ^J\n~dpJk in 
LOCC~ ^-norm). Consequently, one may only hope to make stronger statements about average behaviours. 

This is precisely the general question we address here, being especially interested in the case of high-dimensional 
bipartite quantum systems. We try and quantify in two distinct ways the typical efficiency of the fc-extendibility 
criterion for separability in this asymptotic regime. 

The first approach consists in estimating a specific size parameter (known as the mean width) of the set of 
fc-extendible states when the dimension of the underlying Hilbert space goes to infinity. Comparing the obtained 
value with the known asymptotic estimate for the mean width of the set of separable states then tells us how 
the sizes of these two sets of states scale with one another. The computation is carried out in Section (where 
all needed notions related to high-dimensional convex geometry are properly defined as well) and ends with the 
concluding Theorem |2.5| some technical parts being relegated to Appendix]^ In Sectionthe result is commented 
and comparisons are made between the mean-widths of, on the one hand, fc-extendible states, and on the other, 
separable or PPT states. Besides, a smaller upper bound is derived, in Section and its companion technical 
Appendix on the mean width of the set of fc-extendible states whose extension is required to be PPT (precise 
definitions and motivations to look at this set of states appear there). 

The second approach consists in looking at random mixed states which are obtained by partial tracing over an 
ancilla space a uniformly distributed pure state, and characterizing when these are, with overwhelming probability as 
the dimension of the system grows, fc-extendible or not. Again, comparing the obtained result with the known one for 
separability provides some information on how powerful the fc-extendibility test is to detect entanglement. Section 
introduces all required material regarding the considered model of random-induced states and one possible way 
of detecting their non-fc-extendibility. The adopted strategy is next seen through in Section relying on technical 
statements put in Appendix]^ and concludes as Theorem |6.4| The determined environment dimension below which 
random-induced states are with high probability violating the fc-extendibility criterion is then compared, in Section 
with the previously established ones for violating other separability criteria, and for actually not being separable. 
Finally, generalizations to the unbalanced case are stated in Section(Theorems |8.1| and 8.21, while Section]^ 
exposes miscellaneous concluding remarks and loose ends. 

The reader may have a look at Table for a sample corollary of this study. 

Appendix ^ gathers a bunch of standard definitions and facts about the combinatorics of permutations and 
partitions which are necessary for our purposes. All employed notation on that matter are also introduced there. 
In Appendix the connection is made between computing moments of GUE or Wishart matrices and counting 
permutations having a certain genus. These general observations play a key role in the moments’ derivations of 
Appendices [D| [E| and [Fj which are, as for them, specifically the ones that we need to obtain our various statements. 
To get tractable expressions, though, a formula relating the number of cycles in some specific permutations is 
additionally required, whose proof is detailed in Appendix Appendix finally, is devoted to establishing the 
last crucial ingredient in most of our reasonings, namely bounding the number of non-geodesic permutations (in 
terms of the number of geodesic ones) in some particular instances which are of interest to us. Aside, Appendix 
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Table 1 . Comparison of the average and typical case performance of the fc-extendibility criterion 
with that of the PPT and realignment criteria 


^from the point of view of 

fc-extendibility “beats” ——— 

mean width of 
the set 

entanglement 
detection of 
random states 

PPT 

for fc ^ 11 

for A: ^ 17 

realignment 

7 

for k ^ 5 


[H] is dedicated to proving more precise results than the ones which are strictly needed on the convergence of the 
studied random matrix ensembles, and in generalizing the developed method to establish the asymptotic freeness 
of certain gaussian random matrices. 

Notation. For any Hilbert space H = C”, we shall denote by ’H(H) = 'H(n) the set of Hermitian operators on H, 
and by = 'H+{n) the subset of positive operators on H. For each p G N, we define || • ||p as the Schatten 

p-norm on ■H(H), i.e. || • ||p = (Tr [| • and || • ||oo = limp_,.+oo || • Up- Particular instances of interest are the 

trace class norm || • ||i, the Hilbert-Schmidt norm |j • || 2 , and the operator norm || • ||oo. We shall also denote by 
I?(H) = 'D(n) the set of states on H (positive and trace 1 operators). 

We will in fact mostly consider the case where H = A(g)B = is a (balanced) bipartite Hilbert space. And 

we introduce the additional notations S{A : B) = S{d x d) for the set of separable states on H, P(A : B) = V{d x d) 
for the set of PPT states on H (in both cases in the cut A : B), and for each fc € N, 5fc(A : B) = £k{d x d) for the 
set of fc-extendible states on H (in the cut A : B and with respect to B). 

Preliminary technical lemma. It will be essential for us in the sequel to express in a more tractable way the 
quantity sup^g£^(A-B) Tr(M(T), for any given Hermitian M on A 0 B. Such amenable expression is provided by 
Lemma 11.31 below. 

Lemma 1.3. Let A: G N. For any Mab G 'H{A 0 B), we have 

sup Tr (Mabo-ab) = 

(TAB 

where for each 1 ^ j ^ k, denoting by Idgfc the identity on Bi 0 * * • 0 ^j-i ^ ^j+i 0 * * * 0 we defined 
-^ABfc (j) = ATaBj 0 Idgfc . 

3 

Before proving Lemma [l.3[ let us introduce once and for all the following notation, which we shall later use on 
several occasions: for any Mab^ G FA(A 0 B®^), we define its symmetrisation with respect to B®^ as 

SymA:B'=(^AB'‘) = ^ Yi (Wa 0 f7(7r)B0^AB''(ldA O C/(7r)B0^ 

7TG&{k) 

where for each permutation tt G &{k), U{Tr)Bk denotes the associated permutation unitary on B®^ (see e.g. [5T] for 
further details). 

Proof. By definition, the condition ctab G £fc(A:B) is equivalent to the condition ctab = Tcb^-i SymA.B*" (o'ab'') for 
some (Tab's G V{A 0 B®^). Hence, for any Mab G 'H{A ® B), we have 

sup TrAB [f^ABiTAB] = sup Wab [Mab TrBfc-i SymA-Bf. ((Tab's )] 

<TABe£fc(A:B) o-^gfcGX'(A(giB®'s) 

= sup Ttab's [(Mab ® Wb^-i) SymA.B's((T ab's)] 

''■ab's 6® (Area's) 

= sup Ttab's [SymA:B's (ff^AB 0 IdB's-i)(TAB's] 

'’■ab'sGX'IA^B®'') 

= ||SymA:B's {Mab O Mb's-!) ||^ . 


k 

T slab's (j) 


T = 1 
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Now, for each tt S &{k), (Wa ^ U{Tr)Qk) [Mab 0 Idefc-i) (Wa <8) J7(7r)Bfc)^ = Mab„(i) (8)Idj: 
together, for each 1 ^ j ^ fc, the permutations tt G ©(fc) such that 7r(l) = j, we get 

k 


t(1) 


. Therefore, grouping 


SymA^B*" (^AB ® IdB»=-i) = r X! ^aBj < 8) Id; 


Bf ’ 


i=i 


and hence the advertised result. 


□ 


2. Mean width of the set of Ic-extendible states for “small” k 

2.1. Preliminaries on convex geometry. Let us introduce a few notions coming from classical convex geometry 
which we shall need in the sequel. For any K C H(C") and any M G 'H(C"’) having unit Hilbert-Schmidt norm, 
we define the width of K in the direction M as 

w{K,M)= supTr(MA). 

AeK 

The mean width of K is then defined as the average of w{K, •) over the whole Hilbert-Schmidt unit sphere 
of 7^(0") (equipped with the Haar measure a) i.e. 


sup Tr(MA) 
■*) Laga 


da{M). 


w{K) = I w{K,M)da{M)= f 

7MGSjfs(C") JM&Shs{C” 

This average width w is an interesting size parameter, on its own, but also because it is related to other important 
geometric quantities, such as e.g. the volume radius vrad, which is defined as the radius of the Euclidean ball 
having same volume (i.e. Lebesgue measure). For instance, we have for any convex body K the Urysohn inequality 
w{K) ^ vrad(A'), and for most of the convex bodies K we shall be considering a “reverse” Urysohn inequality 
w{K) ^ /rvrad(iF) for some /r ^ 1. These connections and the precise formulation of these convex geometry results 
are exemplified in Section 

In order to compute the quantity w{K), it is often convenient to re-express it as a Gaussian rather than spherical 
averaging. We thus denote by GUE{n) the Gaussian Unitary Ensemble on C", which is the standard Gaussian 
vector in ?7(C") (equivalently, G ~ GUE(n) if G = {El -\- H^)/y/2 with H a, n X n matrix having independent 
complex normal entries). And we define the Gaussian mean width of iF as 


Wg{K) = EiG^GUEin) 


sup Tr(GA) 
AgA 


Just observing that for G ^ GUE{n), G/||G ||2 is uniformly distributed over and G/||G|| 2 , ||G ||2 are 

independent random variables, we get that the link between both quantities is, setting 7 (n) = EiQ^QUE^n) I|G||2, 
which is known to satisfy 7 (n) n (see e.g. [T], Ghapter 2, for a proof), 

(1) w{K) = -^wg{K). 

l{n) 

Remark 2.1. All the sets K that we will consider in the sequel will actually be subsets o/I?(C"), hence living in 
the hyperplane of EL{C^) composed of trace 1 elements, i.e. in a space of real dimension — 1, rather than n^. It 
would thus seem more natural to define their mean width w{K) as an average width over an^ — 2, rather than — 1, 
dimensional Euclidean unit sphere. The Gaussian mean width wg{K), on the other hand, is an intrinsic notion 
that does not depend on the ambient dimension (because marginals of standard Gaussian vectors are themselves 
standard Gaussian vectors). As a consequence, we see from equation Q that computing the mean width of K as 
if it was a dimensional set is asymptotically equivalent to computing it taking into account that it is in fact a 
— 1 dimensional set. We may therefore serenely forget about this issue. 


Our aim is now to estimate, for any fixed /c G N, the mean width of the set of fc-extendible states on Ac 
A = B = and d —>■ -l-oo. By the definitions above, we have 


) B when 


w{£k{A:B)) = 


1 


-Er 


sup Tr(GABCTAB) 

(Tab (A:B) 


7(^2) 

Using the result of Lemma [l.3| and the notation introduced there, we thus get, 

k 


( 2 ) 


ri;(£fe(A:B)) = 


1 


1{<P) 


E, 


GAB~G(7A(d2) 
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2.2. An operator-norm estimate. As justified above, to obtain the mean width of the set of /c-extendible states 


on 0 C^, what we need is to compute the average operator-norm E 
0 We will show that the following asymptotic estimate holds. 
Proposition 2.2. Fix /c S N. Then, 

k 

EGAB~G;7E(d2) 




for G a GUE matrix on 


y^G*ABfc(j) 

i=i 


d—>-+oo 


2'/kd. 


As a preliminary step towards estimating the sup-norm E 

\2p' 


E-=iGab<=(j) 


, we will look at the 2p-order 


moments E Tr 


G'AB''(j)) 


, p S N, and show that they can be expressed in terms of the 2p-order moments 


of a centered semicircular distribution of appropriate parameter. 

So let us recall first a few required definitions. For any ct > 0, we shall denote by fj-sc'ia^) the centered semicircular 
distribution of variance parameter a^, whose density is given by 


^FSC(a^){x) = ^;^\/4cr2 - X^ll_2a^2cr]{x)dx. 


We shall also denote, for each p G N, by its p-order moment, i.e. It is 

well-known that 

Vp G N, = 0 and = a^^Catp, 

where Catp is the p*^ Catalan number defined in Lemma 


t(p) _ r+°°. 


A.4 


Proposition 2.3. Fix fc G N. Then, when d —>■ -l-oo, the random matrix Gab''(j)^ /d converges in moments 

towards a centered semicircular distribution of parameter k. Equivalently, this means that, for any p G N, 

2p-ij r / \ 2p' 


'E‘GAB~GUE(d^) Tr 


^Gab'=0) 
0=1 


— 0 and Egab~GUE((; 2 ) Tr 


y^GAB^O) 
0 = 1 


d^+aa 


Remark 2.4. Stronger convergence results than the one established in Proposition \2.3\ may in fact be proved, as 
discussed in Appendix\i^ 

_ r /_i. ~ , , \ 

Proof of Proposition\2.3\ Let p G N. Computing the value of the 2p-order moment ETr 


{n'j=i ^AB'^O’)) 


may 


be done using the Gaussian Wick formula (see Lemma B.l for the statement and Appendix B.l for a succinct 
summary of how to derive moments of GUE matrices from it). In our case, what we get by the computations 
carried out in Appendix [P] and summarized in Proposition D.l is that, for any d G N, denoting by tt(-) the number 
of cycles in a permutation. 




7 " ~ \ 




(3) 

ETr 

I X] G'ab''(j) 

= E 

ETr 

]^Gab^(/(^)) 



_Vi=i J 

/:[ 2 p]->.[fc] 



(4) 



= E 

E 

1 

^ 1 


/:[2p]->[fc]AG<p(2)(2p) 

where we defined on {!,..., 2 p}, fp*^^^( 2 p) as the set of pair partitions, 7 = ( 2 p ... 1 ) as the canonical full cycle, 
and for each / : [ 2 p] —>■ [fc], 7 / = 7/=i • • • 7 /=fc as the product of the canonical full cycles on each of the level sets 
of /. 

We now have to understand which A G fp^^^(2p) and / : [2p] —>■ [k] contribute to the dominating term in the 
moment expansion Q, i.e. are such that the quantity tt( 7 “^A) -I- + k — \ Im(/)| is maximal. 

First of all, for any A G fp^^^(2p), we have 

(5) j)(A) -kS( 7 "^A) = 4p- (|A| -f | 7 "^A|) if4p- = 4p - (2p- 1) =2p + l, 

where the first equality is by Lemma A.l| while the second inequality is by equation 


in Lemma A.5 and is an 

equality if and only if the pair-partition A is non-crossing. Next, for any A G fp*^^^(2p) and / : [2p] —)■ [k], we have 


( 6 ) 


+ tt( 7 /= 4p- (^\X\ + 17/< 4p- | 7 ^^| = 4p- (2p- |Im(/)|) = 2p-F |Im(/)|, 
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where the first equality is again by Lemma A.l while the second inequality is by equation (|30| in Lemma A.5 and 


is an equality if and only if the pair-partition A is non-crossing and is finer than the partition of {1,, 2p} induced 
by jf (i.e. / takes the same value on elements belonging to the same pair-block of A). 

Putting equations ([^ and ^ together, we get that for any A S {2p) and / : [2p] —)■ [A:] (just keeping in mind 
that necessarily jj(A) = p), 

(7) tt(7-'A) + tt(77'A) + fc - I Im(/)| ^ 2p + fc + 1, 

with equality if and only if A € NC^'^^{2p) and / o A = /. Since it is well-known that there are Catp elements in 
NC^‘^\2p), and for each of these there are functions which are constant on each of its p pair-blocks, we indeed 
get the asymptotic estimate announced in Proposition |2.3| namely 


ETr 



2p 


AH'" 


(j) 


d—>-+oo 


kP Catp 


□ 


Proof of Proposition^^ The convergence in moments stated in Proposition |2.3| implies that, asymptotically, the 
matrix G'ab'" (j) j /d has a smallest eigenvalue and a largest eigenvalue which are, on average, at most the 

lower-edge and at least the upper-edge of the support of /isc(fe)j i-®- —^Vk and 2^/k. Indeed, convergence of 
all moments of the empirical spectral distribution of ^ab'* (j)) /d entails convergence of all polynomial 

functions, and therefore of all continuous functions with bounded support, when integrated against it. And this in 
turn entails (when applied to continuous functions with support strictly included in [—2y/k, 2'/k]) that the extreme 
eigenvalues of ^ab'=(j)) /d cannot be, on average, strictly bigger than —2\/k or smaller 2\/k. The reader 

is referred to [T], Chapter 2, for all the technical details of the argument. Hence in other words. Proposition 2.3 
guarantees that there exist positive constants Cd -^d^+co 1 such that 


( 8 ) 


E 


y^G*ABfc(i) 


i=i 


^ Cd 2'/kd. 


In the opposite direction. Proposition 


2.3 


only guarantees that the matrix ^ab*" (j)^ /d asymptotically 

has, on average, no strictly positive fraction of eigenvalues strictly below —2Vk or above 2v^. So to show that the 
reverse inequality to Q holds too, a little more care is required. Indeed, to say it roughly, we have to make sure 
that in the moment’s expression the permutations contributing to the non-dominating terms (in d) are not too 
numerous. 

For d G N fixed, it holds thanks to Jensen’s inequality and monotonicity of Schatten norms that 


(9) 


V p G N, E 


Gab*" 


(j) 


< 


E 


g'ab''(j) 


i=i 


2p' 


l/2p 


g; 


ETr 



l/2p 


So let us fix d G N and p G N, and rewrite Q explicitly as an expansion in powers of d, keeping in the sum the 
permutations not saturating equation Q. Being cautious only with the permutations not saturating equation ([^, 
and not with those not saturating equation we get 


( 10 ) 


ETr 


2p 


y^ G^Bk{j) 
0=1 


sG 


L(p+fc)/2j 

E E d-^^ I d^P+^+\ 

yf:[2p]^[k] 


(5=0 


where we defined, for each / : [2p] —>■ [fc] and each 0 ^ d ^ [(p -|- A:)/2J, 

q5g(2p) = {AG'p(2)(2p) : S(77iA)=p+|Im(/)|-2d}. 


In words, tp^^(2p) is nothing else than the set of permutations which have a defect 2d from lying on the geodesics 
between the identity and the product of the canonical full cycles on each of the level sets of /. This justifies in 


particular a posteriori why the summation in (101 is only over even defects (see the parity argument in Lemma 

"Ol. 
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Now, by Lemma G.3 we know that, if 0 ^ i5 ^ L-P/2J) then 

{(/,A) : AG<pg(2p)} 


^ kP Catp X ( ¥ 


26 


And if \p/2] ^ S ^ [{p-^ then trivially 


{(/,A) : XG^f}{2p)} 


^ k'^P 


M! 

2 Pp \ 


^ kP Cat„ X ( ^ 


2 \ P 


Putting everything together, we therefore get, 

2p1 


ETr 


gabi^u) 


vf=l 

1 ^ 


^ kP Catr, 




5=1 



_ ( i ^2p+k+l 


Yet, max |(fcp^/2(i)^‘^ : 1^5<|’p/2]| is attained for <5 = 1, provided p ^ i^dlk^G. So if such is the case. 


V ( — 

^ \ 2d 


5=1 


k 

+ 2 


kp'^^ ^ 


^ p + k k^p'^ ^ k^p^ 


2 Ad? Ad? 


where the last inequality holds as long as p^ k. And hence, under all the previous assumptions, 

2p1 


ETr 


yi gab^o) 




J^(2p) A ^ 2 p + k+l 


So set pd = {2d/k)? for some 0 < e < 1 (which is indeed smaller than {^IdjkYG and bigger than k for d big 
enough, in particular bigger than /c^/^/2). And using inequality in the special case p = pd, we eventually get 


( 11 ) 


E 


Gab*" (j) 


t=i 


< M 


r(2pd) 

^SC(fc) 


1 + 




Ad? ) ) 


^/2pd 


^i+(k+i)/2pd _ 2Vkd. 

d—^-\-oo 


Combining the lower bound in equation ([^ and the upper bound in equation (111 yields Proposition 2.2 □ 


2.3. Conclusion. Combining Proposition 2.2 with equation ([^, we straightforwardly obtain the estimate we were 
looking for, which is stated in Theorem |2.5| below. 


Theorem 2.5. Let fc G N. The mean width of the set of k-extendible states on satisfies 


w{£kiC^:C<^)) 


2 1 

d —^+oo d 


3. Discussion and comparison with the mean width of the set of PPT states 

It was shown in |S] t hat t he mean width of the set of separable states on C‘^ 0 is of order 1 /d^G, And we 
just showed in Theorem 2.5 that, for A: G N fixed, the mean width of the set of /c-extendible states on ® is 
of order 1/d, so that, for d large, 

w{S{C‘^:C^)) < w{£kiC^:C‘^)). 

This result is not surprising: it just means that, when d grows, if k does not grow in some way too, then the set of 
/c-extendible states becomes an increasingly poor approximation of the set of separable states on 0 C'^. There 
had been several evidences, already, in that direction, with examples of highly-extendible, though entangled, states 
(see e.g. [S] and [5S]). 

It is well-known that the exact same feature is actually exhibited by the set of PPT states on 0 whose 
mean width is of order 1/d too. Let us be more precise. 


Proposition 3.1. There exist positive constants Cd, Cd 

e-i/2 


on C‘^ ® satisfies 


1 such that the mean width of the set of PPT states 


Cd 
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Proof. Proposition 3.1 was basically established in [3, but not stated in this exact way and with these exact 
constants, so we briefly recall the argument here for the sake of completeness. 

To get the asymptotic upper bound, we just use 

w(P(C‘^:C^)) < w{V{C‘^ (g) C"^)) 


The last equivalence is a consequence of Wigner’s semicircle law (see e.g. [T], Chapter 2, for a proof) from which it 
follows that 


wg{'D{C^ ® C'^)) — 'E.q^gue{<P) sup Tr(GtT) — 'E,q^que{cP) ||G||oo ~ 2d. 

a-e'D(C‘‘0C‘‘) d-f+oo 

To get the asymptotic lower-bound, we will make use of two results from classical convex geometry. Before 
stating them, we need one more definition: For any convex body K, we denote by vrad(iF) its volume radius, which 
is defined as the radius of the Euclidean ball having the same volume (i.e. Lebesgue measure) as K. 

• Urysohn inequality (see e.g. [32], Corollary 1.4): For any convex body K, 


(12) w{K) ^ vrad(iF). 

• Milman-Pajor inequality (see |26| . Corollary 3): For any convex bodies K,L having the same center of gravity. 


(13) 


vrad(itr n L) vrad(iF — L) ^ vrad(iF) vrad(L). 


Combining equations and ([H, we get that if K, L are convex bodies having the same center of gravity, then 


vrad(iF C L) ^ 


vrad(iF) vrad(L) 
vrad(iF — L) 


vrad(df) vrad(L) 
w{K — L) 


In our case, denoting by F the partial transposition, we have V{d x d) = 
both having the maximally mixed state Id/d^ as center of gravity. Hence, 


vrad(iF) vrad(L) 

^ w{K) + w{L) 

'D{df ) nl?(d^)'", with T>{d^) and 'D{d'^)^ 


, ^ vrad (l?(d^)) 

w{V{d X d)) ^ vrad (iP(d x d)) ^ 2w{V{df)) ’ 

the first inequality being by the Urysohn inequality, and the second being by the Milman-Pajor inequality, after notic¬ 
ing that vrad {V{dfY) = vrad {V{df)) and w{V{dfy') = w{V{d^)). Now, we just argued that w{'D{df)) ~d->.+oo 
2/d, while it was shown in [33] that vrad (V{d^)) '^d^+oo e~^^'^jd. Therefore, 


w(V{d X d)) 


vrad {V{d^)Y 
2w{V{d‘^)) 


e-i/2 

d —>^-|-oo d 


□ 


As a straightforward consequence of Theorem |2.5| and Propositionwe have, roughly speaking, that for /c ^ 11, 
the set of fc-extendible states becomes asymptotically a “better” approximation of the set of separable states than 
the set of PPT states, on average. Indeed, if fc ^ 11, then 2/'/k < so that for d large enough 


4. Adding the PPT constraint on the extension 


The hierarchy of SDPs originally pro posed in m to detect entanglement was in fact slightly different from the 

Indeed, for a given bipartite state pab, the fc**' test would here consist 


1.2 


one that would be derived from Theorem 
in looking for a symmetric extension Pab'" of Pab, while in m it was additionally imposed that this extension had 
to be PPT in any cut of the A: -|- 1 subsystems. This of course increased quite considerably the size of the SDP to 
be solved at each step, but with the hope that it would at the same time decrease dramatically the number of steps 
an entangled state would pass. 

Another hierarchy of SDPs was later proposed in [28] and [29], built on the exact same ideas as those in m- 
It was noticed there that only demanding that the (Bose) symmetric extension of the state be PPT in one fixed 
(even) cut of the fc -|- 1 subsystems already implied a noticeable speed-up in the convergence of the algorithm. It 
therefore seems worth taking a closer look at the set of states arising from these constraints. The latter is properly 
defined as follows. 


Definition 4.1. Let fc G N. A state pab on a bipartite Hilbert space A(8)B k-PPT-extendible with respect to B if 
there exists a state Pab^ on A(8)B®*^ which is PPT in the cut A®B® ^ invariant under any permutation 

of the B subsystems and such that pab = Tr^t-i Pab'"- WA denote by £/^^^(A:B) the set of k-PPT-extendible states 
on A (g) B (in the cut A:B and with respect to B/. 










fc-EXTENDIBILITY OF HIGH-DIMENSIONAL BIPARTITE QUANTUM STATES 


9 


Theorem 4.2. Let fc S N. There exist positive constants Cd 
k-PPT-extendible states on satisfies 


1 such that the mean width of the set of 


Proof. Using the notation introduced in Lemma |1.3| we start from the simple observation that, for any Mab G 
■H(A(g)B), 


sup Tr [Mabctab] = 


sup 


Tr [SymA:B'= {Mab <8 IdBfc-i)crAB'‘] 


^ min (^11 SymA:B'= (^ab 0 Wb^-i) ||^, || [SymA:B'= (^ab <8) IdB''-i)] ^ IL) > 


where F stands here for the partial transposition over the \k/2] last B subsystems, so that in fact 


/ Lfc/2J 


[SymA:Bfc (^AB 8) IdBfc-i)]^ = "I: ( ^ ^aBj 8) Idgfc + ^ Mabj 8) Id 


Bf > 


\ ^'=1 


i=Lfe/2j-Ll 


where F now stands for the partial transposition over B. 

The upper bound in Theorem |4.2| will thus be a direct consequence of the sup-norm estimate 


E, 


GAB'^GUEid^) 


i=i 


d—>-+oo 


V^Vkd. 


The latter is proved in the exact same way as Proposition |2.2| i.e. by first showing that for any p S N, 

2p1 


(14) 


EGAB~G;7£;(<i2) Tr 


u'=i 




and second arguing that also E 


EUGab>^{jY 


J-+C50 J-+00 


ETr 


(E.tlGABdjy 


2p 


l/2p 


. This 


last step will be omitted here since the argument is very similar to the one appearing in the proof of Proposition 
2.2 Concerning the moment estimate (141, it is first of all proved in Appendixthat 


E Tr 


2p' 


G'AB''(j)'^ 


d—>-+oo 


E 


E 


fjih A)+#( 7 j: U)-|-fc-| Im(/)| 


Ae<pG)(2p)/:[2p]^[Lfc/2Jlor[rfc/2l] 


And by the same arguments as the in the proof of Proposition |2.3| we can then identify which A and / actually 
contribute to the dominant order in the latter expression, yielding 


E Tr 


2p 


E ^AB'=(jy 


/ 


d—>-+oo 


E 


\ 


^ fi^P+lk/^j+l ffP+lk/2] 


- 1-1 


/:[2p]^[Lfc/2j] 
\ /o^=/ 


/:[2p]^[rfc/2l] 

foX = f 


d—>+oo 


AGNC(2)(2p) 

Catp ([fc/2jPd2p+L'=/2J+i + (/fc/ 2 ]Pd 2 p+rfe/ 2 Hi^ ^ 


which is the announced moment estimate (141. 


□ 


Comparing Theorem 4.2 to Theorem 2.5 we see that the asymptotic mean width of the set of /c-PPT-extendible 
states is at least \/2 smaller than the asymptotic mean width of the set of fc-extendible states. For instance, the set 
of 2-PPT-extendible states is, on average, asymptotically smaller than the set of 4-extendible states. This however 
does not really shed light on why adding the constraint, at each step in the sequence of tests, that the symmetric 
extension is PPT across one fixed (even) cut would make the entanglement detection notably faster. 

























10 


CECILIA LANCIEN 


5. Preliminaries on random-induced states and witnesses 


We will employ the notation p ~ p,n,s to mean that p = Trc^ with |i/i) a random Haar-distributed pure 

state on (i.e. p describes an n-dimensional system which is obtained by partial-tracing over an s-dimensional 

ancilla space a uniformly distributed pure state on the global “system+ancilla” space). An equivalent mathematical 
characterization of such random state model is p = W/TrW with W ^ Wn,s an (n., s)-Wishart matrix, i.e. W = GG^ 
with G & n X s matrix having independent complex normal entries (see e.g. |34|1. 

Let K C be a convex body. For any p G a standard way of showing that p ^ K is to produce a 

“not belonging to K witness”, i.e. some M G H+(C") which is such that 

supTr(Mcr) < Tr(Mp). 

creK 

By testing p itself as possible such “not belonging to K witness”, we have 
(15) supTr(p(T) < Tr(p^) ^ p ^ K. 

<7GK 

Crucially for the applications we have in mind, the functions p i—> Tr(p^) and p i—> sup^.^^ Tr(pcr) both have nice 
concentration properties around their average. More precisely, we have the two following results. 


Proposition 5.1. Let n, s G N. Then, there exist universal constants c, c' > 0 such that, for any p > 0, first of all 

(|Tr(p') - [Tr(T2)] | ^ 17) < 6"“ + 

and second of all, for any convex body K C I?(C"'), 


P~Mr. 


supTr(pcr) — Et 

aGK 


supTr(Tcr 


^ I? < e 


^—c n sr] 


Proof. To show Proposition 5.1 we will make essential use of a local version of Levy’s Lemma, namely (see |7], 
Lemma 3.4, for a proof): Let Lt C be a subset of the Euclidean unit sphere of R™ satisfying P(n) ^ 7/8. 

Let also / : S'^~^ —)• R be a function whose restriction to Ll is L-Lipschitz and M be a central value for / 
(i.e. P({/ ^ M}) ^ 1/4 and P({/ ^ M}) ^ 1/4). Then, for any 77 > 0, 

P ({1/ - M| ^ 77}) < P (S'™-! \n)+ 
where Cq > 0 is a universal constant. 

It is well-known (see e.g. [31] for a proof) that p ~ pn,s is equivalent to p = XX'^ with X uniformly distributed over 
the Hilbert-Schmidt unit sphere of n x s complex matrices, and the latter can be identified with the real Euclidean 
unit sphere Therefore, one may apply Levy’s lemma above with 17 = {X G : ||A||oo ^ S/i/?!}, 

which is such that P(S^"'®-! \ 17) ^ 6 ““ for some universal constant c > 0 (see e.g. |3], Lemma 6 and Appendix B, 
for a proof). 

Consider first f : X G g'^ns-i Tr((AAl)^), which is 36/n-Lipschitz on 17. Indeed, for any X,Y G 17, 

9 9 

|/(A)-/(y)K (xx^) -{YY^) ^ ^ (||AA^|ioo +llEE^Iloo) (11^112 +rib) 11^-E||2< -||x-y||2. 

The second inequality is just by Holder’s inequality (more specifically ||Ai?G|ji ^ IIAniRlbllGIb) and the triangle 
inequality, after noticing that (XX^)^ - (YY^Y = XX^A + AYY^ with A = X{X^ - pf) + {X - Y)Y^. And the 
third inequality is because, by assumption, for any Z G Ll, nb = 1 and rzl||oo = niL ^ 9/’^- 

Now, the fact that P(r"'®-! \ 17) ^ 6 ““, combined with the fact that |/| is bounded by 1 on implies 

that the average of |/| on \ is bounded by 6 ““, which tends to 0 when s tends to infinity. While the 

Lipschitz estimate for f on LI implies that the average of / on 17 differs from its median by at most C, 
which also tends to 0 when n, s tend to infinity. We can therefore conclude that the average of / is a central value 
of / for n, s big enough. Hence, taking M = E / as central value for /, we get the concentration estimate 


A 


Tr ((Xr)') - EyTr ((Pr)') | > 77 ) 


€ e 


Take next f : X G S'^’"®"! 

|/(A)-/(y)K 


I—>■ supjrg^ Tr(AAf(T), which is G/y^-Lipschitz on 17. Indeed, for any A, F G 17 

6 


sup Tr((AA^ - YY^) a) 

(T^K 


< llAAt _ yytll ^ (1,^11^ + ||y||^) |jy _ y|| 




l|A-y|b. 


The second inequality is just by duality, since K is contained in the unit ball for the 1-norm. The third inequality is 
by the triangle inequality, after noticing that XX^ — EPf = (A — FjAf -|- Y{X^ — Pf). And the fourth inequality 
is by the norm inequality |j • ||oo ^ || • |b because, by assumption, for any Z G 17, ||Z||oo ^ 3f^/ri. 
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Arguing as before, we see that the average of / is a central value of / for n, s big enough (this time, the average 
of I/I on \ is bounded by while the average of / on differs from its median by at most 

Hence, taking M = E / as central value for /, we get the concentration estimate 


A 


supTr(AAV) — Ey supTr(yyV) 

o-Gif o-GA" 


^ V 


^ e-“ + e- 


Hence, we indeed have the two announced deviation probability bounds. 


□ 


K C I?(C") be a convex body. Then, for any rj > 0, 


Combining the two statements in Proposition 5.1 together with equation (15), we get as a consequence: Let 


(16) [Tr(p2)] - E, 


sup Tr (per) 


> V 




where c > 0 is a universal constant. 

From now on, we will in fact consider random-induced states on the bipartite space 0 C^. So let K C 


I?(C“ (g) C“) (such as e.g. 'P(C“:C“) or ffe(C“:C“), k € N). It follows from equation (16) that, for any rj > 0 




(17) E 

where c > 0 is a universal constant. 


supTr(pcr) 

o-ga 


> V 




(p^K)^ ^_g-csmi„(l.dS^)^ 


6. Non fc-EXTENDIBILITY OF RANDOM-INDUCED STATES FOR “SMALL” k 

6.1. Strategy. Our goal in the sequel will be to identify a range of environment size s for which random-induced 
states on ® are, with high-probability, not A:-extendible. In view of equation (0, this may be done by 
characterizing 


s e N : E 




sup Tr(p(T) 

aeSkiC'i-.C'i) 


< E 


'P~Pd2, 


[Tr(p2)] 


Yet by Lemma |l.3| and using the notation introduced there, we have that for any state pab on A ® B, 


sup Tr(pABO-AB) = 

c7AB^^fc(A:B) 


'^Pab'^U) 

i=i 


6.2. An operator-norm estimate. As explained above, to know when random-induced states on A O B are not 

fc-extendible, what we need first is to compute the average operator-norm E X]i=i ^ab'=(j) > for W a {(P,s)- 

. 1—1 . 

Wishart matrix. We will proceed in a very similar way to what was done in Section ^ and establish what can be 
seen as the analogues of Propositions |2.2| and |2.3|but for Wishart instead of CUE matrices. 


Proposition 6.1. Fix A: S N and c > 0. Then, 


E 


1Uab~W^2,^j2 


y^HABfc(j) 

i=i 


As a preliminary step towards estimating the sup-norm E 

\ Pi 


d—>-+oo 


(•\/cfc + 


ELlf^AB^(j) 


, we will look at the p-order moments 


ETr 




, p S N, and show that they can be expressed in terms of the p-order moments of a 
Marcenko-Pastur distribution of appropriate parameter. 

So let us recall first a few required definitions. For any A > 0, we shall denote by Pmp(a) th® Marcenko-Pastur 
distribution of parameter A, whose density is given by 

n I 1 _ j if A > 1 

Pmp{\) X I A/A(a:)da; if A < 1 

where, setting A± = (a/AA 1)^, we defined the function /a by 
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We shall also denote, for each p G N, by its p-order moment, i.e. It is 

well-known that 

p 

r(p) - 

p 

m—l 


tip) 


VpGN, A-Nar; 


where Nar™ is the (p, m)*^ Narayana number defined in Lemma A.4 In particular, = Catp, the Catalan 


number defined in Lemma | A. 41 as well. 

Proposition 6.2. Fix fc G N and c > 0. Then, when d —>■ -l-oo, the random matrix converges 

in moments towards a Marcenko-Pastur distribution of parameter ck. Equivalently, this means that, for any p G N, 


E 


Wab~w^2 ,,^2 Tr 


E^ab'=(j) I 








Remark 6.3. Stronger convergence results than the one established in Proposition \ 6. Sj may in fact be proved, as 
discussed in Appendix\i^ 

Proof of Proposition \ 6. Let p G N. Computing the value of the p-order moment ETr ^ab^O)^ m^iy 

be done using the Gaussian Wick formula (see Lemma B.l for the statement and Appendix B.2 for a succinct 
summary of how to derive moments of Wishart matrices from it). In our case, we get by the computations carried 
out in Appendix [ f] and summarized in Proposition F.l that, for any d, s G N, denoting by !)(•) the number of cycles 
in a permutation. 


ErvAB~Wrf2, Tr 


E^ab4j) 

u'=i 




Tr 


n^AB4/(*)) 


2=1 

-1 


= X! Ewab~W^2„ 
/IpI-s-W 

/:[p]->.[fc]aGe(p) 


where we defined on {1,... ,p}, &{p) as the set of permutations, 7 = (p ... 1) as the canonical full cycle, and for 
each / : [p] —t [A:], 7 / = 7 /=i • • • 7 /=fc as the product of the canonical full cycles on each of the level sets of /. 
Hence, in the case where s = cdf, for some constant c > 0, we have 

\ Pi 


(18) 


EiVAB~w^2,^i2 Tr 


E^ab4j) 

0=1 




= ^ ^ ^ti(«)^2tt(a) + tl(7-'a)+tl(77'a)+fe-|Im(/)|^ 

/:[p]->.[fc] aG6(p) 


We now have to understand which a G S(p) and / : [p] —t [A:] contribute to the dominating term in the moment 
expansion (181, i.e. are such that the quantity 2j)(a) -I- tt( 7 “^a) -I- jl( 7 ^^a) -I- A; — | Im(/)| is maximal. 

First of all, for any a G ©(p), we have 


(19) 


t)(a)-kt)(7 a) = 2p-(|a|-f I 7 a|) ^ 2p - I 7 I = p-t t)(7) = P + 1, 


where the first equality is by Lemma A.l whereas the second inequality is by equation (291 in Lemma A.5 and is 
an equality if and only if a G NC{p). Next, for any a G ©(p) and / : [p] —> [A;], we have 

( 20 ) D(a) -k j)( 7 y ^a) = 2 p - (|a| -f Iq/^a]) ^ 2 p- I 7 /I =p-t j)( 7 /) =p-h |Im(/)|, 

where the first equality is once more by Lemma A.l[ whereas the second inequality is by equation (301 in Lemma 


A.5 and is an equality if and only if a G NC(p) and f o a = f. So equations (191 and (20) together yield that, for 


any a G ©(p) and / : [p] —t [A:], 

(21) 2jl(a) -t tl(7“^a) + ^ “ I Im(/)l ^^p + k+l, 

with equality if and only if a G NC{p) and f o a = f. 


We thus get the asymptotic estimate 



' / k \ 

p- 


( 


ETr 

E^ab^(^’) 


d —^+00 

E E 



. > 


aeAC(p) f-.[p]^[k] 

\ foc^f 

/ 


^2p+fc+l 
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Yet, a function / satisfying / o a = / is fully characterized by its value on each of the tt(a) cycles of a. So there are 
/jtt(a) gucii functions. Hence in the end, the asymptotic estimate 

p-\ 


ETr 


WAB''(j) 


u'=i 


d—>-+oo 


yaeNC{p) 


the last equality being because, for any A > 0, J2aeNC{p) = J2m=i A™Nar™ = 


□ 


6.1\ The argument will follow the exact same lines as the one used to derive Proposition |2.2 


Proof of Proposition 
from Proposition |2.3 

As pointed out there, s how ing the inequality in Proposition |6 .1 1 is easy. Indeed, the convergence in moments 

established in Proposition 


6.2 


implies that, asymptotically, the matrix l^AB'=(j)) /d'^ has a largest eigenvalue 

which is, on average, at least the upper-edge of the support of fiMP{ck)j he. (v^-l-l)^. In other words, it guarantees 
that there exist positive constants Cd —>-d-f+oo 1 such that 


( 22 ) 


E 


wab^u) 

i=i 


^ Cd {Vck + 


Let us now turn to the more tricky part, which is showing the inequality in Proposition 6.1 For d € N 
fixed, it holds thanks to Jensen’s inequality and monotonicity of Schatten norms that 

p \ i/p / r / . \ Pn \ i/p 


(23) 


V p G N, E 


wab*=(j) 

i=i 


E 


y^ hFAB'‘(j) 

i=i 


sG ETr 


y^ hFAB'“0) 
Vi=i 


So let us fix d G N and p G N, and rewrite ( 181 explicitly as an expansion in powers of d, keeping in the sum the 
permutations not saturating equation (21). Being cautious only regarding the permutations not saturating equation 
(201, and not regarding those not saturating equation (191, we thus get the upper bound 

p 


(24) 


ETr 


E^ab4j) 


)]«(. 


L(p+fc)/2j p \ 

E E E d2^’+'=+\ 

/:[p]->[fe] <5=0 m=l J 

where we defined, for each / : [p] —)■ [/c], each 1 ^ m ^ p and each 0 ^ <5 ^ [(p + fc)/2j, 

©/, 5 .m(p) = G 6(p) : t|(a) = TO and jKa) J- (((y^^a) = p J- | Im(/)| - 2d|. 

©/, 5 ,m(p) is thus nothing else than the set of permutations which are composed of m cycles and have a defect 2S 
from lying on the geodesics between the identity and the product of the canonical full cycles on each of the level 
sets of /. This justifies in particular a posteriori why the summation in (241 is only over even defects (see the 
parity argument in Lemma A.2|. Note that the definition of ©/, 5 ,m(p) can actually be extended to all to G N, with 
\&f,s.m{p)\ = 0 if TO ^ p J- I Im(/)| — 2(5, which we shall do in what follows for writing convenience. 

Now, by Lemma G.4 we know that, if 0 ^ <5 ^ Lp/2J j then for any 1 ^ to ^ p. 


|{(/,a) : a G ef,s,mip)}\ < ^E "Nar™ x (2fcV, 


2S 


And if \p/2] ^ S ^ [{p + k)/2\, then trivially for any 1 ^ m ^ p, 


|{(/,a) : aG6/,5,™(p)}| ^ Vp!< K] fc™-^Nar™-d x (2fcV)^ 




What is more, for a given 0^5^ have, making the change of summation index m i—)■ m — e, 


26 


26 




m-eAT^ m-e 


"Narf 


= Y' 








^ (1 + (cfc)™Nar; 

m—l 

= (l + c)2^Mj^Efc)- 
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Putting everything together, we therefore get, 


ETr 


p-| 


E^ab'=(j) 




< 

^ ^^^MP(ck) 


^ ^ 2(1 + c)fcV ^^'^ fc ^ 2(1 + c)feV y I ^2p+fe+i 


Yet, max I (2(1 + : 1 ^ <5 ^ |"p/2] | is attained for (5 = 1, provided p ^ (d/2(l + So if such is 

the case. 


E 

5=1 


+ 2 


2(1 + c)fc^p^\fc /2(1 + c)fc^p^\^ p + A; 4(1 + 4(1 + c)^fc^p^ 


€ 


d2 




d2 


where the last inequality holds as long as p ^ k. And hence, under all the previous assumptions, 

p-| 


ETr 


E U) 

^1 = 1 




So set pd = (d/2(l + for some 0 < e < 1 (which is indeed smaller than (d/2(l + and bigger 

than k for d big enough, in particular bigger than 2{\ + c)k'^^^). And using inequality ( [2^ in the special case p = pd, 
we eventually get 


(25) 


E 


E ^AB''(j) 
1=1 


^ M 


tiPd) 

'■MP{ck) 


(l+ 4(l+^ \y^"'^2+(fc+l)/p. ^ (V^+ 1)2^2. 

\ dd J J d^+oo 


Combining the lower bound in equation (221 and the upper bound in equation (251 yields Proposition 6.1 □ 


6.3. Conclusion. Having at hand the operator-norm estimate from Proposition |6.1| we can now easily answer our 
initial question. It is the content of Theorem |6.4| below. 

Theorem 6.4. Let A: G N, and for any 0 < e < 1/2 define cfik) = (1 — e). Then, there exists a constant 

Ck,e > 0 such that 

[p i £fe(C":C'')) ^ 1 - 

One can take Ck,e = Ce^/k for some universal constant C > 0. 

Proof. As a direct consequence of Proposition |6.1| we have 


E 


Wab~W^2 


sup Tr (Wabctab) 


{Vck + Ifi 2 


d —^-j-CXD k 


d\ 


And since E^ 4 A^yy ^2 TrlT ^d,s^+oo <Ps (see e.g. [T3] or Appendix B.2|, the result we eventually come to after 
renormalizing by TrlT is 


(26) 


E 


P"" 1*^2 =2 


sup Tr(pcr) 

drGffc(C‘i:C‘i) 


1 {yfck + 1)2 2 _ {'/ck + 1)^ 1 


d-f+oo d?cd? 


^d^ = 


ck d? 


On the other hand, Ew~w ^2 ^d,s->.+oo d^s^ + {d'^ fis (see e.g. [13] or Appendix |B.2l), so we also have 

d^{cd^fi + {d^ficd^ 


Ep'-P,2„,2 [Tr(p )] [d^cdffi 




Now, if c = (1 — e){k — 1)2/4A: for some 0 < e < 1/2, then 

(V^-fl)^ _ (^1+ E = 

ck 


4e 


< 


8e 


(/c-l)(l-e) A:-l' 


So by equation (171, we have in such case 

{p i £fc(C‘^:C^)) ^ 1 - exp [-Ckd^d^{tlkd^f) = 1 - exp {-Cd^t^jk) , 


^ P~Pd2,ad-^ 

for some universal constant C > 0. 


□ 
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7. Discussion and comparison with other separability criteria 


For each A: G N, define Ck-ext as the smallest constant c such that a random state p on ® induced by an 

environment of dimension c(P is not fc-extendible with high probability when d is large. That is, 


^k—ext — illf 1 C 




d —^-j-cxu 


What we established in Theorem 6.4 is that Ck-ext ^ {k — 1)^/Ak. 

Yet, we know from |7] that for c > 0 fixed, p ~ PcP^ctP is with high probability entangled when d —> +oo: the 
threshold for p ^ PdP^s{d) being with high probability either entangled or separable occurs for some s{d) = so{d) 
with d^ < so{d) < d^ log^ d. So what we proved is that if c < (/c — l)^/4fc, i.e. if A: > 2c + 2a/c(c + 1) + 1, then this 
generic entanglement will be generically detected by the A:-extendibility test. 

Furthermore, it is well-known (see e.g. |34| 1 that p ~ PdP^d? is equivalent to p being uniformly distributed on 
the set of mixed states on ® (for the Haar measure induced by the Hilbert-Schmidt distance). As just 
mentioned, when d —>■ -|-c», such states are typically not separable. Now, for A; ^ 6, {k — 1)^/4A: > 1, so such 
states are also typically not A:-extendible. Hence, entanglement of uniformly distributed mixed states on 
is typically detected by the A:-extendibility test for k ^ 6. 

Let us define, in a similar way to what was done for the A:-extendibility criterion, Cppt, resp. Cra, as the smallest 
constant c such that a random state p on ® induced by an environment of dimension is, with probability 
tending to one when d tends to infinity, not satisfying the PPT, resp. realignment, criterion. We know from [?] 
that Cppt = 4, whereas we know from |1] that Cra = (8/37r) . Now, for k ^ 17, {k — 1)^/4A; > 4, and for k ^ 5, 
{k — 1)^/4A: > (S/Stt)^. So roughly speaking, this means that the A:-extendibility criterion for separability becomes 
“better” than the PPT one at most for k ^ 17, and “better” than the realignment one at most for A; ^ 5. This is to 
be taken in the following sense: if k ^ 17, resp. A: ^ 5, then there is a range of environment dimensions for which 
random-induced states have a generic entanglement which is generically detected by the A:-extendibility test but not 
detected by the PPT, resp. realignment, test. 

Note also that for the reduction criterion [22], it was established in [^ that the threshold for a random-induced 
state on either passing or failing it with high probability occurs at an environment dimension d, hence 

much smaller than for all previously mentioned criteria. 

8. The unbalanced case 

For the sake of simplicity, we previously focussed on the case where H = A 0 B is a balanced bipartite Hilbert 
space. One may now wonder what happens, more generally, when A = and B = with cJa and de being 
possibly different. It is easy to see that the results from Theorems |6.4| and |2.5| straightforwardly generalize to the 
case where dA and de both tend to infinity (but possibly at different rates). The corresponding statements appear 
in Theorem 18.II below. 


Theorem 8.1. Let A: G N and let dAjde G N. The mean width of the set of k-extendible states on 
(with respeet to ) satisfies 

dA.dB—>-+00 'fkvdj^d'Q 

Also, when dA, dn —> -koo, a random state on 0 whieh is sampled from p,dAdB,cdAdB> with c < {k — 1)^/4A:, 
is with high probability not k-extendible (with respect to ). 

Oppositely, when one of the two subsystems has a fixed dimension and the other one only has an increasing 
dimension, the sets of A:-extendible states with respect to either the smaller or the bigger subsystem exhibit different 
size scalings. This is made precise in Theorem |8.2| below. 


Theorem 8.2. Let A: G N and let dA,dB G N. If dA is fixed, the mean width of the set of k-extendible states on 
0 (with respect to ) satisfies 

dB^-kOO yLf \fdAdfi 

Whereas if dB is fixed, the mean width of the set of k-extendible states on (with respect to ) satisfies 

2C{dB,k) 1 




>+oo y/k \/ dAds 
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Proof. Using the same notation as in the proof of Proposit ion |2.3[ we start in both cases from the exact expression 
for the 2p-order moment (slightly generalizing Proposition D.ll 

\ 2 p‘ 

"Hi: 

U'=l / /:[2p]->[fc]AG‘p(=)(2p) 


(27) 


EGAB-'G;7E(<iAdB) Tr 


First, fix d^- The argument then follows the exact same lines as in the proof of Proposition |2.3| Indeed, the pair 
partitions A £ *p(^)(2p) contributing to the dominant order in ds in the expansion (271 are the Catp non-crossing 
pair partitions A S NC^‘^\2p), for which j)( 7 “^A) = p-l- 1. Moreover, for each of these A, the functions / : [2p] —)■ [fc] 
contributing to the dominant order in ds in the expansion (27) are the functions which are such that / o A = /, 
for which tt( 7 j^A) -I- A: — | Im(/)| = p + k. So we eventually get 


EGAB~G;7E(dAdB) Tr 



dB—>- + 00 


Catp 


Now, fix de • Again, the pair partitions A S (2p) contributing to the dominant order in d\ in the expansion 
(27) are the Catp non-crossing pair partitions A € (2p), for which tt( 7 “^A) = p -I- 1. So consider one of these 

A. Observe that, for any 0 ^ d ^ L(P + ^)/2jj if / : [2p] —>■ [fc] is such that there are exactly S pair blocks of A on 
which / takes two values, then necessarily (((y^^A) -\- k — \ Im(/)| ^ p -|- fc — 2d. Indeed, the case d = 0 is already 
known. So let us describe precisely what happens in the case d = 1, i.e. when there is exactly 1 pair block of A on 
which / takes 2 values. 

• If amongst these 2 values, at least 1 of them is also taken on another pair block of A, then there exist transpositions 
r, t' and a function g satisfying g o \ = g, such that 7 / = ^qTt' and | Im(/)| = | Im(p)|. Hence, 

I 1 ( 77 ^A) = tt(TV 77 ^A) ^ tt( 77 ^A) - 2 = p -f I Im(p)| - 2=p+\ Im(/)| - 2. 


• If none of these 2 values is also taken on another pair block of A, then there exist a transposition r and a function 

g satisfying go \ = g, such that 7 / = 7 gr and | Im(/)| = | Im(( 7 )| -|- 1. Hence, 

tt( 77 ^A) = tl(r 77 U) ^ 11 ( 77 ^A) - 1 =p-k |Im(p)| - 1 =p-h |Im(/)| - 2 . 

And this generalizes in a similar way to d > 1. Yet, for a given 0 ^ d ^ [(p-|-fc)/2j, there are (^)fc^(fc— functions 

which take 2 values on exactly d pair blocks of A (assuming of course that p ^ k). So we eventually get 


EGAB~GC/E(dAdB) Tr 


G'AB'=(d) j 

U=i / 


2p 


/L(p+fe)/ 2 j 

^ Catp F ^ ( ^ ) (fc - I d^+'^dP+i 

\ d =0 


/L(p+fe)/ 2 j 






= Catp kP [1 + 


\ s=o 
k-1 


d| 


L(p+fc)/ 2 j 




One can then argue as in the proof of the derivation of Proposition |2.2| from Proposition |2 . 3| that we additionally 

. . l/ 2 p 

k 


have E 


~ limp_>+oo (ETr Gab^ j j > when either dn -t +00 or dA 

automatically yields the two announced statements on the mean width of 
Remark 8.3. In the situation where du is fixed, if we had an exact expression 

V p e N, E(3^g..,GC/E(dAdB) Tr GABfc(j)J . Af(dA,p), 

then we would be able to conclude without any further argument that 

k 


- 1 - 00 . This 

□ 




E 


gab^ 

i=i 


^ lim M{dx,pY^‘^^. 

dA-H-oo p-H-00 
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This would indeed follow from the convergence result of |20| for non-commutative polynomials in multi-variables 
with matrix coefficients (in our case, variables and dg x dg coefficients). 


The asymmetry in the definition of fc-extendibility appears more strikingly in this unbalanced setting. Indeed, 
for a finite fc S N, a given state on A 0 B may be fc-extendible with respect to B but not fc-extendible with respect 
to A. It is only in the limit k —)■ +oo that there is equivalence between the two notions: a state on A 0 B is 
A:-extendible with respect to B for all fc G N if and only if it is fc-extendible with respect to A for all A: G N (and if 
and only if it is separable). 

However, what Theorem |8.1| stipulates is that, even for a finite A: G N, when both subsystems grow, being k- 
extendible with respect to either one or the other are two constraints which are, on average, equivalently restricting. 
On the contrary, what Theorem |8.2| shows is that when only one subsystem grows, and the other remains of fixed 
size, being A:-extendible with respect to the bigger one is, on average, a tougher constraint than being A:-extendible 
with respect to the smaller one (as one would have probably expected). 

This is to be put in perspective with some of the original observations made in m- It was indeed noticed that 
checking whether a state on is A:-extendible with respect to requires space resources which scale as 


d' -\-k — l 
k 


)d 


when implemented. It was therefore advised that in the unbalanced situation of d “big” and d' “small”, 

one should check A:-extendibility with respect to rather than C^, the former being much more economical. 
On the other hand, it comes out from our study that, in this case, an entangled state is likely to fail passing the 
A:-extendibility test for a smaller k when the extension is searched with respect to than when it is searched with 
respect to . But understanding the precise trade-off seems out of reach at the moment. 


9. Miscellaneous questions 

9.1. What about the mean width of the set of A:-extendible states for “big” k? All the statements proven 
sofar, regarding either the A:-extendibility of random-induced states or the mean width of the set of A:-extendible 
states, converge towards the same (expected) conclusion: for any given A: G N, the A:-extendibility criterion becomes 
a very weak necessary condition for separability when the dimension of the considered bipartite system increases. 
So the natural question at that point is: what can be said about the A:-extendibility criterion on ® when 
k = k{d) is allowed to grow in some way with d? Unfortunately, most of the results we established rely at some 
point on the assumption that A; is a fixed parameter, and therefore do not seem to be directly generalizable to the 
case where k depends on d. 

There is at least one estimate though that remains valid in this setting, which is the lower bound on the mean 
width of A:-extendible states. 


Theorem 9.1. There exist positive constants C 4 -^d ^+00 1 such that, for any k{d) G N, the mean width of the set 
of k{d)-extendible states on 0“^® satisfies 

u; (C":C‘^)) 

Proof. Let d, k{d) G N. For any p G N, the exact expression for the 2p-order moment established in Proposition 
|D.l| of course remains true. So by the same arguments as in the proof of Proposition |2. 3 1 we still have in that case 
at least the lower bound 

2 p- 

^ Y rf«(7-'A)-Ltt(77A)-hfc(d)-|Im(/)| 

/:[2p]^[fe(d)]AG'p(2)(2p) 

^ Catp A:(d)Pd2p+'=(‘^)+^ 

This lower bound on moments in turn guarantees, as explained in the derivation of Proposition |2.2| from Proposition 

□ 1 such that we have the inequality 

k{d) 

(j) 

i=i 


EGAB~GC/E(d2)Tr 


/ k{d) ^ \ ■ 

I (j) 1 


2.3 that there exist positive constants Cd 


E, 


GAB~GUE{d^) 


^ Cd 2yjk{d)d, 

which yields the announced lower bound for the mean width of £k{di){C‘^'.C^). 


□ 


Theorem 9.1 only provides a lower bound on the asymptotic mean width of 5fc(C“:C“) when k is allowed to 
depend on d. It is nevertheless already an interesting piece of information. Indeed, as mentioned in Sectio n [3] 
we know from that the mean width of the set of separable states on is of order 1/d^/^. Theorem 


9.1 
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therefore asserts that, on 0 C^, one has to go at least to k of order d to obtain a set of /c-extendible states whose 
mean width scales as the one of the set of separable states. 

Furthermore, it may be worth mentioning that the proof of Proposition |2.2| actually provides additional informa¬ 
tion, namely an upper bound on the mean width of /c-extendible states which remains valid for a quite wide range 
of k. 

Theorem 9.2. For any d, k{d) £ N, provided that k{d) < and d is big enough, the mean width of the set of 
k{d)-extendible states on satisfies 

f k{dy/^ Ind' 




VHd)d 


exp 


I d2/5 


Proof. Let d, k{d) £ N with k{d) < . Taking pd = (dlkipfiffil^ in equation (which is indeed, as required, 

bigger than k{d) and smaller than (2d/A;(d))^/^ for d big enough) we get 


E, 


GAB~GC/E(d2) 


k{d) 

^ G'AB'=(j) 
i=l 


< 2^/fc (^d exp ind + 


2d2/5 


In d -1- In - 
4 


The latter quantity is smaller than 2yjk{d) d exp (fc(d)^/^ In d/d^/®) for d big enough, which yields the advertised 
upper bound for the mean width of □ 


Of course, the upper bound provided by Theorem 9.2 is interesting only for k{d) < ko{d) = d^/^/(ln d)5/7. 


Nevertheless, since the set of fc-extendible states contains the set of /c'-extendible states for all k' k, we also have 
as a (potentially weak) consequence of Theorem |9.2| that for fc(d) ^ fco(d) and d big enough, 

2e (In d)®/^'^ 


(4w(C‘':C")) <u;(£fc„(,)(C'':C'')) < 


d8/7 


Theorems 9.1 and 9.2 together imply in particular the following: in the regime where k{d) grows with d slower 
than d itself, both the ratio w;(iS(C‘^:C‘^))/u>(£fc(d)(C^:C‘^)) and the ratio w(£k(d)iC‘^'C‘'‘))/ w(T>{C^ ® C^)) are 
unbounded. To rephrase it, for fc(d) having this growth rate, the set of fc(d)-extendible states lies “strictly in 
between” the set of separable states and the set of all states from an asymptotic size point of view. 


9.2. When is a random-induced state with high probability fc-extendible? The result provided by Theorem 


by partial tracing on C® a uniformly distributed pure state on C“(8)C“®C® is with high probability not fc-extendible. 
But what can be said about the case s > s{k,d)? Or more generally, can one find a reasonable s'{k,d) such that 
if s > s'(fc, d), then a random mixed state on ® obtained by partial tracing on C® a uniformly distributed 
pure state on 0 C® is with high probability fc-extendible? 

By the arguments discussed in extensive depth in [7], one can assert at least that there exists a universal constant 
c > 0 such that ckdf log^ d is a possible value for such s'(fc, d). We will not repeat the whole reasoning here, but let 
us still give the key ideas underlying it. 

Define £fe(C^:C'^) the translation of £fc(C‘^:C^) by its center of mass, the maximally mixed state Id/df, i.e. 

£k{C'^-.Cy = i^p-^^ : p££:fc(C'':C^)|. 


6.4 is only one-sided: it tells us that if s < s{k^ d) = dr[k — 1)^/4A:, then a random mixed state on C“ 0 obtained 


Define also i5fc(C'^:C‘^)° the convex body polar to £fc(C‘^:C'^), i.e. 

£k{C‘^:Cy° = {A £ n{d^) : V A £ £fc(C‘^:C^), Tr (AA) < 1} . 
What then has to be specifically determined is (see [7], Section 2, for further comments) 


£ N : E 




sup Tr 

AG£fc(C‘i:C‘i)° 


Id , 

^ d2 I 


< 1 


One can first of all use the fact that, roughly speaking, when d, s —>■ -l-oo, the random matrix p — ld/df for p ^ Pd'^.s 
“looks the same as” the random matrix Gjdf^/s for G ~ GUE^d?) (see |7], Proposition 3.1 and Remark 3.2 as 
well as Appendices A and B, for precise majorization statements and proofs). In particular, there exists a constant 
C > 0 such that, for all d, s £ N with (say) d^ ^ s ^ d^, we have the upper bound 

, _ sup Tr ( ( p - ^ I A I < -^-J= E.a^auEid'^) _ sup Tr (GA). 

AG£fc(C‘i:C‘i)° \\ ^ J J U V® AGffclC'iiC'i)” 
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Next, due to the fact that, again putting it vaguely, the convex body is “sufficiently well-balanced” (see 

|7], Section 4 as well as Appendices C and D, for a complete exposition of the f'-position argument), we know that 
there exists a constant C" > 0 such that, for all d € N, we have the upper bound 


E, 


G~GUE(d^) 


Ag£ 


sup 


Tr (GA) 



GUE{<P) sup Tr (GA) < G'd'^ log d. 

AeSkiC<‘:C‘‘) J 


Now, EG^( 3 ( 7 £;((i 2 ) sup^g£^((^d.(^d) Tr (GA) is nothing else than the Gaussian mean width of £’fc(C‘^:C^), which is 

namely wg (£’fc(C'^:C‘^)) '-^d^+x 2d/y/k. 

Putting everything together, we see that 


the same as the Gaussian mean width of £fc(C“:C“), so for which we have an estimate thanks to Theorem 


2.5 


AG£fc(C‘i:C‘i)» 



^G 


y/kd log d 


for some constant G > 0 independent of d,s,k £ N, which implies as claimed that if s > G^/cd^log^d, then 
supAe£,(c£c<i)° Tr ((p - Id/d^) A) < 1. 


Remark 9.3. Let us briefly comment on a notable difference, from a convex geometry point of view, between the k- 
extendibility criterion and other common separability criteria. In the case of k-extendibility, computing the support 
function ofS^ is easier than computing the support function of its polar while for other separability relaxations 
it is usually the opposite. Indeed, for a given traceless unit Hilbert-Schmidt norm Hermitian A on (g) C^, we 
have for instance the closed formulas 


'-£fc(C£C-i) 


(A) = sup{Tr(AA) : A £ £’fe(C'^:C‘^)} = 


A 
IIA^ 


/ip(c.:c<i)o(A) =sup{Tr(AA) : A £ P(C‘^:C'^)°} = d^ 
whereas the dual quantities (^) cannot be written in such a simple way. 

This explains why in the case of Sk it is the mean width that can be exactly computed, contrary to the threshold 
value which can only be approximated, while for other approximations of S the reverse generally happens. 
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Appendix A. Combinatorics of permutations and partitions; short summary of standard facts 

Let p £ N. We denote by &{p) the set of permutations on {!,... ,p}. For any tt £ G{p), we denote by jl(7r) 
the number of cycles in the decomposition of tt into a product of disjoint cycles, and by |7r| the minimal number 
of transpositions in the decomposition of tt into a product of transpositions. We also define 7 £ &{p) as the 
canonical full cycle {p . ■. 1). More generally, we shall say that c is the canonical full cycle on a set {A, • • ■, ip} with 
ii <■■■< ip if c = {ip ... ii). 

Some standard results related to &{p) are gathered below (see e.g. [SD], Lectures 9 and 23, for more details). 
Lemma A.l. For any tt £ &{p), tt(7r) -|- |7r| = p. 

Lemma A. 2. d : (tt, c) £ &{p) x &{p) >-->• |7r“^c| defines a distance on ©(p), so that for any 7r,ij £ 6(p), 

(28) k“^7r| -I- |7r| = d{g, tt) -|- d{id, tt) ^ d(id, ^) = k|, 

with equality in ( |28[ ) if and only if n lies on the geodesic between id and c. And whenever this is not the case, there 
exists S £ {1,... ,p — 1} such that -I- |7r| = k| -|- 26. 
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Definition A.3. A partition of {1 ,... ,p} is a family A = {Ii,..., Ii,} of disjoint non-empty subsets of {1 ,... ,p} 
whose union is {1,... ,p}. The sets /i,... ,/l are called the blocks of X. If each of them contains exactly 2 elements, 
A is said to be a pair partition of {1,... ,p}. We shall denote by ^{p) the set of partitions of {1,... ,p}, and by 
the set of pair partitions of {l,...,p}. Note that = % if p is odd. Remark also that, whenever 

p is even, the set of pair partitions of {1 ,... ,p} is in bisection with the set of pairings on p} (i.e. the set 

of permutations on p} which are a product of p/2 disjoint transpositions). We shall therefore make no 

distinction between both. 

A partition of {1,... ,p} is said to be non-crossing if there does not exist i < j < k < I in {1,... ,p} such that 
i,k belong to the same block, j,l belong to the same block, and i,j belong to different blocks. We shall denote by 
NC{p) the set of non-crossing partitions of {l,...,p}, and by NC^^'^{p) the set of pair non-crossing partitions of 
{1,... ,p}. Note that NC^‘^\p) =% if p is odd. 

A well-known combinatorial result regarding non-crossing partitions is the following. 


Lemma A.4. The number of non-crossing partitions of {1,... ,p} and the number of pair non-crossing partitions 
of {1,..., 2p} are both equal to the p*^ Catalan number 


Cato = 


p-Il 


More precisely, for any 1 ^ m ^ p, the number of non-crossing partitions of {1,... ,p} which are composed of 
exactly m blocks is equal to the {p, Narayana number 


Nar™ = 

^ p+l 


1 fp + l 


p-l 
m — 1 


Obviously, these numbers are such that J2m=i = Catp. 

With these definitions in mind, we can now state a special case of particular interest of Lemma |A.2[ 
Lemma A.5. Denote by 7 the canonical full cycle on {1,... ,p}. Then, for any tt S &{p), 

(29) | 7 "V|-f | 7 r| ^ I 7 I =_p-1, 


with equality in (291 if and only if tt lies on the geodesic between id and 7. The latter subset of &{p) is in bijection 
with the set of non-crossing partitions of { 1 ,... ,p} (by the mapping which associates to a given partition the product 
of the canonical full cycles on each of its blocks). We shall thus write tt G NC{p) in such case, not distinguishing a 
geodesic permutation from its corresponding non-crossing partition. 

More generally, let {Ii,... ,Ii/\ be a partition of {1,... ,p} and denote by the canonical full cycles on 

Ii , ... ,Il. Then, for any tt G &{p), 

( 30 ) 1(71 • • ■ 7 L)”^ 7 r| -f | 7 r| ^ I71 • •.7^] = p - L, 

with equality in ( 30 ) if and only if tt lies on the geodesic between id and ■■■The latter subset of &{p) is in 
bijection with the set of non-crossing partitions of ... ,p} which are finer than L U • • • U Ir, which itself is in 
bijection with NC{\Ii\) x • • • x NC{\Il\). 

Combining Lemma |A.5 with Lemma A.4| we can in fact say the following: Let G &{p) and assume that 
its decomposition into disjoint cycles is = ci • • • cl where, for each 1 ^ ^ L, Ci is of length pi (hence with 
Pi + ■ ■ ■ + Pl = p)- Then, for any tt G &{p), + kl ^ p — L, and 

|{7re6(p) : |ij“^7r|-I-|7r| = p — L}| = Catpj X • • • X Catp^. 

Having this easy observation in mind might be useful later on. 


Appendix B. Computing moments of Gaussian matrices: Wick formula and genus expansion 

When computing expectations of Gaussian random variables, a useful tool is the Wick formula (see e.g. [3S] or 
m, Lecture 22, for a proof). 

Lemma B.l (Gaussian Wick formula). Let Xi, ..., Xg be jointly Gaussian centered random variables (real or 
complex). 

If q = 2p -\- 1 is odd, then E [Xi ■ ■ ■ Xq] = 0. 

p 

If q = 2p is even, then E [Xi ■ ■ ■ Xq] = 
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B.l. Moments of GUE matrices. A first important application of Lemma [B.l| is to the computation of the 
moments of matrices from the Gaussian Unitary Ensemble. Indeed, for any g G N, we have 

EG~GU£;(n) Tr (C^) = X/ ® [^h.h ■ ■ ■ , 

l^/l . ,lg^n 

where the Gij, 1 ^ i,j ^ n, are centered Gaussian random variables satisfying E[Gi^jGi'j'] = So what 

we get applying the Wick formula is that, for any p € N, 

EG^GUE(n)Tr(G2P+l) =0 

Eg^GUE(„) Tr (G^P) = ^ 

AG*p(2)(2p) 


where for each pair partition A = {{*i, ji}, ■ • •, {*p, jp}} of 2p}, b(A) is the number of free parameters 

h, ■ ■ ■ ihp € {1,...,n} when imposing that V 1 ^ to ^ p, h^+i = Ij^, ^,^+1 — Identifying the pair partition 
{{zi, ji},... ,{zp,ip}} with the pairing (fy ji)... (ipjp) and denoting by 7 the canonical full cycle ( 2 p ... 1), the 
latter condition can be written as V 1 ^ z ^ 2p, fy-iA(i) = h- So in fact, b(A) = tt( 7 “^A) and the expression above 
becomes 

Eg-^GUEW Tr (G^^) = ^ 

AGqjG)(2p) 

We thus have the so-called genus expansion (see e.g. m, Lecture 22) 

Lp/ 2 J 

Eg^gc/eh Tr (G^P) = ^ P(2p, 

5=0 


where for each 0^5^ Lp/ 2J > we defined P(2p, d) as the number of pairings of {!,..., 2p} having genus 6 , i.e. 
P{2p,S)= {Ae^(2)(2p) . |(^-iA)=p+1-2(5} = {A 6^3(2)(2p) : tt(7-iA) + tt(A) = 2p + 1 - 25} . 


Equivalently, P{2p, S) is the number of pairings of {1,..., 2p} having a defect 25 of being on the geodesics between 
id and 7 . Hence, P{2p, 0) is the number of pairings of {!,..., 2p} lying exactly on the geodesics between id and 7 , 
i.e. the number of non-crossing pair partitions of {l,...,2p}. So P(2p,0) = |A^G^^^(2p)| = Gafy, and we recover 
the well-known asymptotic estimate 

EG^GUE(n) Tr (G^P) ^ CafyrzP+i. 


B.2. Moments of Wishart matrices. A second important application of Lemma [B.l| is to the computation of 
the moments of matrices from the Wishart Ensemble. In such case, a graphical way of visualising the Wick formula 
has been developed in m, to which the reader is referred for further details and proofs, a brief summary only being 
provided here. 

In the graphical formalism, a matrix X : C™ —)■ C" is represented by a “box” with two “gates”, one specifying 
the size to at its entrance and the other specifying the size n at its exit. For X : C"* —>■ C" and Y : C" —>■ C”^, the 
product XY : C" —>■ C" is represented by a wire connecting the exit of Y to the entrance of X. For Z : C™ —C"*, 
the trace Tr(Z) is represented by a wire connecting the exit and the entrance of Z. 

Let W be a (rz, s)-Wishart matrix, i.e. W = GG^ with G a nx s matrix with independent complex normal entries. 
Representing by ♦ a n-dimensional gate and by T a s-dimensional gate, the quantity Tr(WP) is then graphically 
represented by p boxes G and p boxes G^ connected by wires in the following way. 



For any a G &(p), we will denote by ffa the diagram obtained from the one above by “erasing” the boxes, just 
keeping their gates, and then connecting, for each 1 ^ i ^ p, the entrance of the z*^ box G to the exit of the a(z)*^ 
box G^, and the exit of the z*^ box G to the entrance of the a(z)*^ box G^. Doing so, loops connecting 

rz-dimensional gates and tl(a) loops connecting s-dimensional gates are obtained. And the graphical version of the 
Wick formula tells us that 


Ti(lVP) = H 

ae&(p) ctGS(p) 
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In the special case where s = n, this can be rewritten as a so-called genus expansion (see e.g. m) 

Lp/2 J 

Ew^w„,„Tt{WP) = Y. S{p.6)nP+^-^\ 
s=o 

where for each 0 ^ J ^ Lp/2Jj we defined S{p,d) as the number of permutations on p} having genus S, 

i.e. ^(p, (5) = |{a e 6(p) : tt(7“^a)-I-tt(a) = p-|-1 - 2i5}|. Since {a G 6 (p) : tl(7“^a)-I-tt(a) = p-|-1} = fVC'(p), 
we have S'(p, 0) = Catp and hence recover the well-known asymptotic estimate 

„ Tr (lE^) ~ Catp . 

’ n—>-+00 


Appendix C. One needed combinatorial fact: relating the number of cycles in some specific 

PERMUTATIONS ON EITHER [p] X [k] OR [p] 


Let a G 6(p) and / : [p] —t [k]. We define d/ on [p] x [fc] as 


V {i,r) G [p] X [k], af{i,r) = 


(a(z),/(a(i))) if r = f{i) 
(z,r) if r 7^ f{i) 


We also define 7 on [p] x [A] as (7, id), where 7 G ©(p) is the canonical full cycle (p ... 1). 

We would like to understand what is the number of cycles in 7“^d/. For that, it will be convenient to do a bit 
of rewriting. Let us first extend the definition of dy and 7 to [p] x ({0} U [/c]). We shall denote by ctf and 7 the 
respective extensions. Note that since / takes values in [fc], we have af{i, 0) = (i, 0) for all i G [p]. 

We will now make two easy observations. 


Fact C.l. For any / : [p] —>■ [k\, define for each i G [p], as the transposition on [p] x ({0} U [k]) which swaps 
(i, 0) and (i, f{i)), and set fif = We then have, for any a G ©(p), 

d/ = where V (z,r) G [p] x ({0} U [A:]), a'{i,r) = 


{a{i),r) if r = 0 
(z,r) ifr 7^ 0 


The advantage of expressing d/ in this way is that a' is particularly simple: it acts as a x id on [p] x {0} 
and does nothing on [p] x [k]. Furthermore, due to the cyclicity of jl(-), a direct consequence of Fact C.l| is that 
jl(7“^d/) = tt(7j ^d'), where 7/ = It may then be easily checked that 7/ decomposes into fc -|- 1 disjoint 

cycles as stated in Fact |C. 2] below. 


Fact C.2. For any / : [p] — >■ [A:], we have 

7/ = Cl • • • CfcC, 

with c = ((p, /(p))... (1, /(I))), and for each r G [A:], Cr = ((p, Sr{p)) ■ • • (1, Sr(l))), where for each i G [p], Sr(i) = 0 
fii) = r and Sr(i) = r if f [i) 7 ^ r. 


Example C.3. For the sake of concreteness, let us have a look at a simple example. In the case where p = 
4, A: = 3, and f is defined by /(I) = /(2) = /(4) = 1 , /(3) = 2, we obtain that the cycles in fif are Ci = 
((4,0)(3,1)(2,0)(1,0)), C2 = ((4,2)(3,0)(2,2)(1,2)), C3 = ((4, 3)(3,3)(2, 3)(1, 3)), c = ((4,1)(3, 2)(2,1)(1,1)). This 
is schematically represented in Figure^ where the elements in Ci are marked by “i”, for i G {1,2,3}, and the 
elements in c are marked by 


Figure 2. / : [4] ^ [3] such that f-\l) = {1,2,4}, f-\2) = {3}, /-i(3) = 0. 


r G {0} U [3] 


3 

3 
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2 
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Lemma C.4. Let / : [p] — ^ [A:] and define 7/ G &ip) as 7/ = 7/=i • • where for each r G [k], ')f=r is the 

canonical full cycle on f~^{r). Then, for any a G &{p), 

(31) + 1 + A: - I Im(/)|. 


Proof. In 7y there are, first of all: 

• A; — I Im(/)| cycles of the form ((l,r)... {p,r)) for r G [A:] \ Im(/), because for any {i,r) G [p] x ([A:] \ Im(/)), 
jj'^a'{i,r) = 7j ^(i,r) = (7-1(1), r). 

• 1 cycle ((!,/(!)) ■ • ■ {p,f{p))), because for any i G [p], y^a'{i, f{i)) = f'J^{i,f{i)) = (7"H0./(7"H*)))- 

For the cycles in which belong to none of these two categories, there are two crucial observations to be made. 

First, for any i,j G [p], {i, 0) and (j, 0) belong to the same cycle of 7j^d' if and only if i and j belong to the same 

cycle of And second, for each i G [p] and r G [A:] \ {/(a(i))}, there exists j G [p] such that (*,r) belongs to 

the same cycle of 7^^d' as (j,0). Indeed, for any i G [p\, we have on the one hand 

7-ia(i) = (7-i)^+ia(i), with L = inf{Z ^ 0 : /((7-i)'+ia(i)) =/(a(z))}. 

While we have on the other hand. 


/(7 ^a(O) =/(«(*)) ^ 7/^a'(z,0) = (7 1 q;(z), 0) = (7^ ^^(z), 0), 

-V,)) . ^ iyo^i^L-1, (77V)'(z,0) = ((7-i)VM*),/(«W)) 

\(77V)^(z,0) = ((7-i)Va(z),0) = (77M*),0) 

So there are in fact exactly remaining cycles in 77^d'. 


□ 


Example C.5. Looking at the same example as before, namely p = A, k = 3, and f such that f~^(l) = {1,2,4}, 
/-i(2) = {3}, /-i(3) = 0, we see that, for a = (14)(23), the cycles in are: 

• ((1, 3)(2,3)(3,3)(4, 3)), because 3 ^lm.{f). 

• ((1,1)(2,1)(3, 2)(4,1)), because /(I) = 1, /(2) = 1, /(3) = 2 and /(4) = 1. 

• (1,0) and ((2, 0)(2, 2)(1,2)(4, 2)(3, 0)(4, 0)(3,1)), corresponding to the cycles (1) and (2,3,4) in^^^a. 

Putting together these preliminary technical results, we straightforwardly obtain Proposition |C.6| below. 


Proposition C.6. Let / : [p] — t [A:] and define 'jf G &{p) as 7^ = 7/=i • • • "ff^kj where for each r G [A;], 7/=r- is the 
canonical full cycle on f~^{r). Then, for any a G S(p), 

(32) j)(7-id/) = 11(7710)-hA:-|Im(/)|. 

Proof. This is a direct consequence of Lemma |C.4[ just noticing that = |1(7-1 q;/) + 1. □ 


Appendix D. Proof of the moments expression for “modified” GUE matrices (mean-width of the 

SET OF A;-EXTENDIBLE STATES) 


The goal of this Appendix is to generalize the methodology described in Appendix |B.1| in order to compute the 
2p-order moments of the matrix X)j=i ^ab'" (j)- Recall that this issue arises when trying to estimate the mean 
width of the set of A:-extendible states. We are thus dealing here, not with standard GUE matrices, but with 
d^-dimensional GUE matrices which are tensorized with d*^-i-dimensional identity matrices. 

For any ii,..., Z2p G [A;], we can write 

X! ^ab*" {ii)LX2 ■ ■ ■ ^ab-" i^2p)r^^Xi > 

U ,...42pe 


Tr 




JJGab'=(u) 


i=i 


where for each j G [2p] and each Ij = (oj, 6 ],..., b^) , Ij+i = (oj-i-i, 5]^^,..., 6j+i) G [d]^^^, we have 

f^AB^ {ij\j ,lj + i ^ (a j,bf + + i * 

Consequently, for any / : [2p] —)■ [A;], we have 

■ 2p 


Tr 


J^Gab^ (/(*)) 


i=l 


ai,...,a2pG[d] bi,...X2peIj 




where If = jfci, ... ,b 2 p G [d]'^ : V z G [ 2 p], V r G [A:] \ {/(z)}, = 6 }|. 
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What we therefore get by the Wick formula for Gaussian matrices is that, for any / : [2p] —>■ [A:], 

■ 2p 


EGAB~G;7E(<i2)Tr 


]^Gab'^(/{ 0) 


= dKA)+KA/)^ 

AG^p(2)(2p) 


where for each pair partition A = ji}i ■ ■ ■ j {ipj jp}} of 2p}, b(A) is the number of free parameters 

ai,...,a2p S [d] when imposing that V 1 ^ m ^ p, a^^+i = and b(A/) is the number of 

free parameters bi,...,b 2 p & If when imposing that V 1 ^ to ^ p, As noticed 

before, identifying the pair partition {{*i, ji}, • ■ •, {ip,jp}} with the pairing {ii ji) ... (ipjp) and denoting by 7 the 
canonical full cycle (2p ... 1), the latter conditions may be written as 


7-1(1; 


V i G [2p], = at and y r G [k] 

So in fact, b(A) = tt(7“^A) and b(A/) = tt(7“^A/), where 

V (i,r) G [2p] X [fc], 7(7r) = (■y(i),r) and Xf{i,r) = 


K-^%)=K ifr = /W 


= if r ^ f{i) 


(A(i),/(A(i))) if r = /(z) 
{i,r) if r ^ f(i) 


What is more, we know by Proposition C.6 that jl(7“^A/) = jl(7^^A) -I- fc — | Im(/)|, where 7/ is the product of the 
canonical full cycles on the level sets of /. 

Let us summarize. 

Proposition D.l. For any d G N and any p G N, we have 



7 ^ ^ 


'7- 

EGAB~GC/E(d2)Tr 

E^AB'=(j) 

= E EGAB~GC/E(d2)Tr 

n^AB^ (/(*)) 


_V=i / 

/:[2p]->-[fc] 



/:[2p]^[fe]AG'P<2)(2p) 

Appendix E. Proof of the moments expression for partial transposition of “modified” GUE 

MATRICES (mean-WIDTH OF THE SET OF A:-PPT-EXTENDIBLE STATES) 

The goal of this Appendix is to compute the 2p-order moments of 'Yf=i G’ab'=(j)^j where P stands here for the 
partial transposition over the [fc/2] last B subsystems. Recall that this issue arises when trying to estimate the 
mean width of the set of fc-PPT-extendible states. 

Using the same notation as in Appendix]^ and reasoning in a completely analogous way, we have that, for any 
/ : [2p] -t [k], 

■ 2p, 


Tr = Y E 

ai,...,a2pG[d] bi,...,b2pelj 

i if f{i) < [fc/2j 
z -I- 1 if /(z) > [A:/2J 

What we therefore get by the Wick formula for Gaussian matrices is that, for any / : [2p] 

J AGq3C)(2p) 




where for each 1 ^ z ^ 2p, Xi = 


and Xi = 


z -I- 1 if /(z) < [A:/2J 

* if fii) > L^/2J 


[k], 


EGAB~G(7E(d2)Tr 


where for each pair partition A = {{Uj ji}, ■ • ■, {*pj jp}} of {!,..., 2p}, b(A/) is the number of free parameters 
61,..., b 2 p G [d]^ when imposing that for all z G [2p], first = 6” if r /(*)> and second the one condition 

if /(0:/(-^(0) ^ L^/2J or /(z),/(A(z)) > [/c/2j, while the two conditions d:^-E(!) = and 
bi\Y^ = if/(*) ^ Lfc/2J, /(A(z)) > Lfc/2J or /(z) > Lfc/2J, /(A(z)) < [k/2\. 

Let us rephrase what we just established. Fix A G tpi^i(2p). For functions / : [2p] —>■ {1,..., [fc/2j} = [[/c/2j] 
or / : [2p] —>■ {[fc/2j + 1,... ,k} = [|’A:/2]], the number of free parameters associated to the pair (A,/) is the 
same as the one observed in Appendix On the contrary, for functions / which are not of this form, extra 
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matching conditions are imposed. So these will for sure not contribute to the dominating term in the expansion of 

,2p-| ^ 

into powers of d. Consequently, we have the asymptotic estimate 


ETr 


(E,tiGAB<=(jr 


'^GAii~GUE{cP) Tr 


2 p 


U'=l 


d—>-+oo 


E 


E 


^#(7 ^>')+tt(7/^A)-|-fc-| Im(/)| 


A6<p(2)(2p) /:[2p]->[Lfc/2j]or[rfc/2]] 


Appendix F. Proof of the moments expression for “modified” Wishart matrices (A:-extendibility 

OF RANDOM-INDUCED STATES) 

The goal of this Appendix is to generalize the methodology described in Appendix |B.2| in order to compute 
the p-order moments of the matrix ^ab*" (j)- Recall that this issue arises when trying to characterize k- 

extendibility of random-induced states. We are thus dealing here, not with standard Wishart matrices, but with 
(d^, s)-Wishart matrices which are tensorized with (i^“^-dimensional identity matrices. 

Representing by • a d-dimensional gate and by T a s-dimensional gate, the matrix Wab'=( 1)) for instance, may 
be graphically represented as in Figure 


Figure 3. Wf^-Qk{l) — Aab'"(1)Aab'‘(1)^j with Aab'"(1) — Gabi C)IdB2...Bfc 

-^AB''(1) -^AB''(1)^ 



The products Wab'" (1) Wab'" (1) and Wab'" (1) WAb*" (2), for instance, are then obtained by the wirings represented 
in Figure]^ 


Figure 4. Wab''(1)Rab'=(1) (on the left) and Wab'= (I)Rab'' (2) (on the right) 



So what we get by the graphical Wick formula for Wishart matrices is that for any f : [p] ^ [k ], 




n^AB4/(*)) 


2 = 1 


ae&(p) ctGS(p) 


7 “/)gtt(a) 


where d/ is defined by 


af : {i,r) £ [p] x [k] i-D- 


(a(i),/(a(i))) if r = f{i) 
{i,r) if r ^ f{i) 


and where 7 stands for 7 applied to the first argument. Indeed, for each a G &{p), there are tl(7“^Q;) loops connecting 
the d-dimensional gates corresponding to A, tl(7“^d/) loops connecting the d-dimensional gates corresponding to 
Ri,... and tl(Q;) loops connecting s-dimensional gates. This is because for each 1 ^ ^ p, on subsystems A 
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and the entrances (respectively the exit) of the box Xj^^k{f(i)) are connected to the exits (respectively 

the entrance) of the a{iY'^ box X^^k{f{a{i))Y. 

What happens in the special case p = 2 and fc = 2 is detailed in Figures]^ andbelow as an illustration. 


Figure 5. /(I) = /(2) = 1. On the left, a = id: 





Finally, we also know by Proposition C.6 that for any a € ©(p) and / : [p] —>■ [fc], denoting by 7y the product of 
the canonical full cycles on the level sets of /, we have |1(7 “^q;/) = tt(7j^a) + fc — | Im(/)|. 


Putting everything together, we eventually come to the result summarized in Proposition |F.1| below. 
Proposition F.l. For any d, s G N and any p G N, we have 


EiTAB-Wrf2_,Tr 


= X! ^Tr 

n^AB'= (/(*)) 


/ J 

f-[p]->-[k] 



— ^ ^ £^#(7 ^a)+t)(7^ V)+fe-| lm(/)|^tt(a)^ 

/:[p]->-[fc]aG6(p) 

Appendix G. Counting geodesics vs non-geodesics pairings and permutations 

Let us r ecall once and for all two notation that we will use repeatedly in this section, and that were introduced 
For any p, to G N with to ^ p, we denote by Catp = the p*^ Catalan number, and by 


A.4 


in Lemma 

Nar™ = (m-i) (P’ to)*^ Narayana number. 

G.l. Number of pairings of 2p elements which are not on the geodesics between the identity and the 
canonical full cycle. 

Lemma G.l. Let p € N and denote by 7 the canonical full cycle on {1, ... ,2p}. For any 0 ^ d ^ \p/2\, define 
the set of pairings having a defect 26 of being on the geodesics between id and 7 as 


(2p) = {Ae^(")(2p) : tt(7-iA)=p + l-2d}. 


>( 2 ) 


Then, the cardinality of^^^\2p) is upper bounded by Catp (p'^/4)'*. 


To prove Lemma G.ll (and later on Lemma G.3l we will need the simple observation below. Roughly speaking, 

* ' /OA 

it will allow us to assume without loss of generality that, in the decomposition of an element of T2p) into p 
disjoint transpositions, the ones “creating” the 26 geodesic defects are the 26 first ones. 
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Fact G.2. Let be a permutation on {1,..., q} and ti, T 2 , T 3 he three disjoint transpositions on {1,, q}, for some 
integer qj? 6. Define = ‘^ti, = gTiT 2 , = gTiT 2 T 3 , and assume that 

(33) = tt(0 + 1, = «(0 + 2, = tt(^) + 1. 

Then, there exists a permutation tt of the three indices {1,2,3} such that, defining this time = <;T7r(i), = 

^ ^ ^ '^7r(l)'^7r(2)'^7r(3) ^ have 

tt(.rW) = +1, = tt(0, = tt(0 +1- 


Proof. Assume that and ti = (ii ji),T2 = { 1232 ),= {hjs) satisfy equation (33l. This means that ii,ji belong 
to the same cycle of i 2 ,j 2 belong to the same cycle of and 13, J3 belong to two different cycles of So let 
us inspect all the scenarios which may occur. 

• Cl cfc\ and C 2 ^*24^^ cfc^, with ci,C 2 two different cycles of <j: If is G cf and js G c\ then the re-ordering 

1, 3,2 is suitable. If is G cf and js G c| then the re-ordering 2, 3,1 is suitable. If is G dj and js G C2, for a,b G {x, y}, 
then both re-orderings 1,3, 2 and 2, 3,1 are suitable. And similarly when the roles of is and js are exchanged. 

• c ^*1-4^^ c'c" *-* 24 ^^ with c a cycle of while c^ = c" and c' ^*24^^ c^c^: If is G (f and js G then the 

re-ordering 2,3,1 is suitable. If is G c“, for a G {x,y}, and js G then the re-ordering 1,3,2 is suitable. And 
similarly when the roles of is and js are exchanged. □ 


As an immediate consequence of Fact G.2 we have the following: Let ? G &(2p) and X = ti- ■ - Tp G *}3^^^(2p). 


Define for each 1 d, q d, p, = g n ■ ■ ■ Tq, as well as = c- Assume next that, for some 0 < J < L(P + tl(‘?))/2j, 

Then, there exists a permutation tt of the p indices {1,... ,p} such that, defining this time for each 1 ^ q ^ p, 
Tn.(q), as well as = ?, we have 

- 1 if g G {2e : 1 < e < 5} 

-I- 1 if g ^ {2e : 1 ^ e ^ 5} 


(q) 


(34) 


V 1 < g < p. 


Since A = ri • • • Tp = T7r(i) • ■ • Tt^{p), we see that we may always assume without loss of generality that, given the 
transpositions ri,..., Tp in the decomposition of A are ordered so that A is under the canonical form (34). The 
behaviour of the function g G [p] j)(c^^^) under this hypothesis, depending on the value of 6, is represented in 

Figure]^ (in the special case p = 6 and {((g) = 1). 


Figure 7. Case p = 6 and tt(<j) = 1. From left to right: (5 = 0, (5=1, i5 = 2 and (5 = 3. 



With this result in mind, let us now turn to the proof of Lemma |G.1| 

Proof of Lemma \G.l\ Given A = {ii ji) ■ ■ ■ (ipjp) G 5p^^^(2p), we will always assume from now that the transposi¬ 
tions (ii ji),..., {ipjp) in its decomposition are ordered so that A is under the canonical form ( [M| for 7“^. This 
means the following: defining, for each 1 ^ g ^ p, the permutation = 7~^(*i Ji) ■ • • (*g Jq) and the integer 
n(g) = j)(A^*^), as well as A^°^ = 7“^ and n(0) = D(A^°^) = 1, we have, for any 0^5^ L7'/2j5 

[n(g) = n{q — 1) — 1 if g G {2e : 1 < e < (5} 

I n(g) = n{q - 1) -|- 1 if g ^ {2e : 1 < e < (5} 


A G ^ (2p) V 1 < g < p. 


In particular, A G NC^^\2p) ^ V 1 ^ g < p, n(g) = n{q — 1) -I- 1, and we know that there are precisely Catp 
possibilities to build such pairing A. This implies that, for each 1 ^ (5 ^ \p/‘2‘\^ there are necessarily less than 
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(^2^) ■ ■ ■ ^ Catp- 2 S possibilities to build a A G Indeed, for the choice of the 25 first disjoint 

transpositions we can use the trivial upper bound that would consist in picking them completely arbitrarily, while 
the p — 25 last ones have to be chosen so that they form a non-crossing pairing of the 2p — A5 not yet selected 
indices. Now, we just have to observe that 

^2{p - 25 
2 


1 ) 


X Catp—25 — 


2p • • • (2p - 4(5-k 1) 


(2p-4(5)! 


1 

4S 


22-5 (p-2(5)!(p-2(5-k 1)! 

p\{p+l)\ {2p)\ 


{p - 25)\{p - 25 + l)\ p!(p-l-l)! 


P 


^ 2'2S ^ Catp, 


which completes the proof. 


□ 


G.2. One needed generalization: bounding the number of pairings of 2p elements which are not on the 
geodesic path between the identity and a product of (few) cycles. The proof of Proposition 2.2 crucially 


relies at some point on a statement of the same kind as the one appearing in Lemma |G.1| Nevertheless, what 
we actually need there is a slight generalization of the latter. More specifically, we have to bound the number of 
pairings which have some defect of lying on the geodesics between the identity and, not only a full cycle, but also a 


product of (few) cycles. So let us give the following extension of Lemma G.l which is really directed towards the 
application that we have in mind. 

Lemma G.3. Let p G N. For any f : [2p] —^ [fc] and any 0 ^ (5 ^ L(P + |I™(/)I) /2J; define the set of pairings 
having a defect 25 of being on the geodesics between id and jf (the product of the canonical full cycles on each of 
the |Im(/)| level sets of f) as 

q3^5(2p) = {AG‘p(2)(2p) : UlfX)=p+\lmif)\-25}. 

Then, for any 0 ^ (5 ^ \p/2\, we have the upper bound 

{(/,A) : AGq3g(2p)}| . 

Proof. We will follow the same strategy and employ the same notation as in the proof of Lemma |G.1| Given 
A = (*i ji) ■ • • {ipjp) G and / : [2p] [k], we will always assume that the transpositions (zi jfi),..., {ipjp) 

in the decomposition of A are ordered so that A is under the canonical form p4| for This means the following: 

defining, for each 1 ^ g < p, A^"^^ = 'yj^iii ji) ■ ■ ■ {iq jq) and n{q) = (KA^*^), as well as A^°^ = and n(0) = 
jj(A(o)) = |Im(/)|, we have, for any 0 < (5 ^ [{p + |Im(/)|)/2j, 


(35) 


A G ^p](2p) V 1 g < p. 


n(g) = n{q — 1) — 1 if g G {2e : 1 ^ e ^ 5} 
n(g) = n{q — 1) -I-1 if g ^ {2e : 1 ^ e ^ 5} 


In particular, A G ip^^Q(2p) V 1 ^ g ^ p, n(g) = n(g — 1) -I- 1, and we know that there are precisely Catp 
possibilities to build a pair (/, A) satisfying this condition (because the latter holds if and only if both constraints 
A G NC^'^\2p) and / o A = / are fulfilled). 

In the case Lp/2Jj notice that after the 25 first steps, we are left with a permutation c having |Im(/)| 

cycles, and we have to impose that the partial pairing A = (*25+1 j25+i) • • • (^2p j2p) lies on the geodesics between id 
and c. Now, the number of such partial pairings is the same as the number of partial pairings lying on the geodesics 
between id and 7j, for any function / : [2p] —)■ [A] whose level sets are the supports of the cycles of g. Hence, to build 

a pair (/, A) meeting our requirements, we have at most x Catp_25 possibilities. Indeed, 

for the 25 first disjoint transpositions we can use the trivial upper bound that would consist in picking them, as well 
as the values of / on them, completely arbitrarily, while for the p — 25 last ones we have to impose that they are 


non-crossing and that / takes only one value on a given transposition. Now, we know from the proof of Lemma G.l 
that (2 P)... (2(p-/+i))catp_25 ^ (p^/4)‘^Gatp. So we get as announced that there are less than Gatp(p^/4)^ 
pairs (/, A) satisfying condition (35l. □ 


G.3. One needed adaptation: bounding the number of permutations of p elements which are not 


on the geodesic path between the identity and a product of (few) cycles. The proof of Proposition 6.1 


requires a statement analogous to the one appearing in Lemma |G.3| but for permutations instead of pairings. In 
order to derive it, we need first to explicit a bit how an element of &{p) can be put in one-to-one correspondence 
with an element of {2p) whose pairs are all composed of one even integer and one odd integer. 
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To a full cycle c= {ii ... Ji) on {1,..., we associate the pairing Ac = (2zi 2z2 — 1) • • • {2ii 2zi — 1) on {1,..., 21}. 
The reverse operation is obtained by collapsing the two elements 2i and 2 j — 1 to a single element i for each 
1 ^ i ^ p. Then as expected, we associate to a general permutation a = ci ■ ■ • Cm £ &{p) the pairing Aq, = 
Aci---Ac„ 

Observe that, denoting by 7 the canonical full cycle either on {1,... ,p} or on 2p}, we have 

(36) V a G 6 (p), tl(a) + tl(7"^a) = S(7"^A„). 

Indeed, the cycles of 7 “^Aq, are precisely cycles of the form {2ii ... 2ii) for (i; ... ii) a cycle of a (supported on even 
integers) and of the form (2b' — 1 ... 2ii — 1) for (b' ... ii) a cycle of 7 “^a (supported on odd integers). So what 
equation ( [M| shows is that the elements of &{p) having a given geodesic defect are in bijection with the elements 
of *p(^)(2p) with even-odd pairs only and having the same geodesic defect (between id and 7 in both cases). In 
particular, we recover the well-known bijection between NC{p) and NC^'^\2p) (because a non-crossing pairing is 
necessarily composed of even-odd pairs only). 

Next, for any function g, either from [p] to [fc] or from [2p\ to [b], we will denote by ■jg the permutation, either on 
{1,... ,p} or on {1,..., 2p}, which is the product of the canonical full cycles on the level sets of g. For any function 
f ■ [p] —>■ [b], we define the function / : [2p] —)■ [fc] by f{2i) = f{2i — 1) = /(z) for each 1 ^ z ^ p. It is then easy to 
see that we have more generally 

V / : [p] ^ [k], V a G 6 (p), ))(«) + 

This simple observation will allow us to derive, as a slight adaptation of Lemma |G.3| a corresponding estimate 
for permutations instead of pairings. 


Lemma G.4. Let p G N. For any / : [p] —> [k], any 0 ^ A ^ [(p -|- | Im(/)|) /2J, and any 1 ^ m ^ p — 26, define 
the set of permutations which are composed ofm disjoint cycles and which have a defect 26 of being on the geodesics 
between id and 7 / (the product of the canonical full cycles on each of the |Im(/)| level sets of f) as 

&f,s,m{p) = {a e &{p) : tl(a) = TO and j)( 7 ^^a)-|-j)(a) = p-|-| Im(/)| - 25}. 

Then, for any 0^5^ \p/'2\ o,nd any 1 ^ to ^ p — 25, we have the upper bound 

25 

|{(/,a) : a€ ©//,m(p)}| ^ (4A:V)'^ 

e =0 


Proof. We just observed that, for any 0^5^ Lp/ 2J and 1 ^ to ^ p — 25, the following equivalence holds 


(37) 


aG6/,5,m(p) ^ tl(a) = TO and Aq, G fp 


( 2 ) 

/A’ 


where fp~(. denotes the set of pairings having a defect 25 of lying on the geodesics between id and 7 ?, as defined in 

r.o J 

Lemma IG.3I 

In particular, a G 6 /p_m(p) tt(a) = to and Aq, G 


and we know that there are precisely fc” 


Nar” 


possibilities to build a pair (/, a) satisfying this condition (because the latter holds if and only if the three constraints 
j)(a) = m, a G NC{p) and f o a = f are fulfilled). 

For the case 1^5^ [p/2j, we will mimic the proof of Lemma G.3 So let (/, a) be such that a G ©/, 5 ,m(p) and 
assume without loss of generality that the transpositions (ii ji ),..., (ipjp) in Aq are ordered so that Aq is under the 
canonical form (34) for 7 j. This means that the partial pairing {i 2 S+i j 25 +i)''' {ipjp) is on the geodesics between 
id and some c with | Im(/)| cycles, and the number of such partial pairings is the same as the number of partial 
pairings being on the geodesics between id and some 77 with | Im(/)| = | Im(/)|. Hence, to count how many ways 
there are of constructing what happens on {zi, ji,... ,i 25 ,j 2 s}, we have the trivial upper bound that would arise 
if picking the 25 first transpositions in Aq, as well as the values of / on them, completely arbitrarily. This yields 
a number of possibilities of at most k"^^ ■ While on {z 25 -i- 1 j j 2 ( 5 +i, ■ ■ ■ fi 2 p,j 2 p}, we have to impose 

that the p — 26 last transpositions in Aq are non-crossing, and that, when collapsed into a permutation of p — 25 
elements, the latter has between to — 25 and to cycles and the function / takes only one value on each of them. 
This leaves us with a number of possibilities of at most Putting everything together, we see 


that the number of pairs (/, a) satisfying condition 

1 ) 


is less than 


ui5 


(2{p-6- 


25 


( 2 p' 


26 

,25 


e =0 


e=0 


which is exactly what we wanted to show. 


□ 
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Remark G.5. The upper bound we established in Lemma G.4 is probably far from optimal (e.g. it is likely that 
the exponent A5 in the polynomial pre-factor in k and p can be improved). But this does not really matter for our 
specific goal. Nonetheless, in the special case of non-geodesic permutations between id and 7 on {1,... ,p}, it is in 
fact quite easy to obtain an upper bound which scales as p^^ for the ratio between the number of 25 non-geodesic 
permutations with a given number of cycles and the number of geodesic permutations with the same number of cycles. 
We present the result in Lemma \G.6\ below, the problem being that the proof method does not seem to generalize so 
straightforwardly to the case that we truly need, that is the one of non-geodesic permutations between id and 7 /. 

Note also that very similar looking upper bounds had previously been derived regarding the cardinality of the set 
@s{p) = {a G S(p) : tt( 7 “^a) + tt(a) =p+ 1 — 2(5}, which is the union of the sets Gg^mip) defined in Lemma G.6, 
for 1 ^ m ^ p. In particular, it was established in m, Lemma 12, that for any 0^5^ Lp/2J, 

| 65 (p)| < \6o{p)\p^^ = Ca,tpp^^. 

However, this is definitely even less enough for our purpose: the latter really requires an upper bound on the number 
of permutations which have a given defect and a given number of cycles in terms of the number of permutations 
which have no defect and the same (or a related) number of cycles. 

On the other hand, one may have hoped for a stronger result than these simply counting ones. For instance 
something like 

(i(id, ( 2 ) + (i(a, 7 ) = d(id, 7 ) + 2(5 ^ 3 a' : d{a,a') = 25' and d{id,a')-\-d{a',p) = d{id,p), 

with S' ^ 95 and with the mapping (p : a i-G- a' satisfying |(/)“^(( 2 ')| ^ p''^, for some coefficients 9 ,k. However, 
determining whether this kind of statement holds or not seems to remain an open question. 

Lemma G.6. Let p G N and denote by 7 the canonical full cycle on For any 0^5^ Ll'/2J o^nd 

1 ^ m ^ p—25, define the set of permutations which are composed of m disjoint cycles and which are 25-away from 
the geodesics between id and 7 as 

&s,mip) = {a G &{p) ■ ‘iia)=m and jl( 7 “^( 2 ) + jl(a) = p + 1 - 2(5}. 

Then, the cardinality of & 5 ^rn{p) is upper hounded in terms of the cardinality o/©o,m(p) os 


I ©(5,m {p) 1^1® 0,m ip)\ 


Proof. Let p S N and 1 ^ m ^ p. We know from m, Theorems 4.1 and 4.2, that there exist polynomials Pq of 
degree q, for 0 ^ q ^ Lp/ 2J ) snch that for any 0 ^ (5 < [(p — m)/2j, 

s 


(38) 


|65.m(p)| = 


p! 


p + 1 - 25 
22-5 (2(5)! I m 


E 

£=0 


P- 1 
TO — 1 + 2e 


P^{m)Ps-e{p + 1 - TO - 25). 


What is more, one can check from the explicit expression provided there for the polynomials Pq that, for any x ^ 0, 
Pq{x) ^ (2x)®. As a particular instance of equation (38), we have 

^p + 1\ / p - 1' 

TO 7 V TO — 1 / 


| 6 o.m(p)| = p! 

And as a consequence, we get by a brutal upper bounding that, for any 1 ^ ^ ^ [(p — to)/2 J, 


|65.m(p)| 

|® 0 .m(p)| 


< 


< 


1 

22-5(25)1 

1 

22-5(25)1 

1 

22-5(25)1 

3 \ -5 


m— 1 

n 

2^0 

E 

e=0 


p + 1 — 25 — Z -p-r . P — ™ — J 


P + 1 — Z 


e =0 j =0 


p — m 


TO 


2 £ 


m-\- j 

2 '' m" 2 ^“*^ (p + 1 — TO — 25 )"^“ 


P,:{m)Ps-e{p + 1 - TO - 25) 


-5 .^3S 


X (5+l)2V 


And this is precisely the claimed upper bound. 




□ 


Remark G.7. There is a close link between the problem we are concerned with and the one of finding tractable 
expressions for the so-called connection coefficients of the symmetric group (the reader is referred e.g. to |18| for more 
on that topic). Glosed formulas are actually known for the connection coefficients o/©(p), involving the characters 
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of its irreducible representations. But unfortunately, they are not really handleable in there full generality. And it 
seems it is only in some specific cases that more manageable forms can been obtained (i.e. in the first place as a 
sum of positive terms, so that one can see more easily what its order of magnitude is). The two situations which 
are well-understood are, on the one hand when the function f is constant (which corresponds to the case where 7 / 
is the canonical full cycle, and hence has a particularly simple cycle type, that is treated e.g. in m), and on the 
other hand when the defect 26 is 0 (which corresponds to the case of so-called top connection coefficients). 

Appendix H. Extra remarks on the convergence of the studied random matrix ensembles 

For any Hermitian M on C", we shall denote by Ai(M),..., \n{M) G R its eigenvalues, and by Nm its eigenvalue 
distribution, i.e. the probability measure on R defined by 


NM = -y^6, 

71 




In words, for any / C R, Nm{I) is the proportion of eigenvalues of M which belong to I. 

H.l. “Modified” Wishart ensemble. Fix /c G N and c > 0. Then, for each d G N, let VF ^ 'AidP,c(P define 
the random positive semidefinite matrix Wd on by 

1 - ^ 

(39) 




i=i 


Proposition 6.2 establishes that when d —>■ + 00 , the eigenvalue distribution of Wd converges in moments towards a 
Marcenko-Pastur distribution of parameter ck. But a stronger result actually holds, namely that there is convergence 
in probability of towards p.MP{ck)- What is meant is made precise in Theorem 


H.l 


below. 


Theorem H.l. For any J G R and any e > 0, 


d—>-+oo 

where the matrix Wd is as defined in equation 


lim ~ TMP(ck){I) \ > s) — 0, 


Theorem 

Proposition 


H.l 


6.2 


is a direct consequence of the estimate on the p-order moments ETr Wab'^ (j) 


from 


combined with the estimate on the p-order variances Var Tr Wab'^ U)^ from Proposition 

H.2 below. The proof, which follows a quite standard procedure, may be found detailed in [T| and sketched in [?]. 


Proposition H.2. Let p G N. For any constant c > 0, 




y^IFABfc(j) 




= o (d 2 p+'=+i) . 

d —>'-|-oo 


Proof. Let p G N. We already know that ^ETr ^ab*" (j) j 

Proposition |6.2| Consequently, the only thing that remains to be shown in order to establish Proposition |H.2| is 

p -| \ 2 / _ •^2 




) thanks to 


that we also have E (Tr 


ELl W^AB‘ 


^ combinatorics involved in the 


proof of the latter estimate is very similar to the one already appearing in the proof of the former. We will therefore 
skip some of the details here. 

To begin with, let us fix a few additional notation. We define 71 = (p ... 1) and 72 = (2p ... p + 1) as the 
canonical full cycles on {1,... ,p} and {p + 1,..., 2p} respectively. Also, for each functions fi : {1,... ,p} —)■ [fc], 
/2 : {p + 1 ,..., 2 p} [k], we define the function /i ,2 : [ 2 p] -t [fc ] by /i ,2 = fi on { 1 ,... ,p} and /i ,2 = /2 on 

{p + 1,..., 2p}. Then, by a slight generalization of Proposition C. 6 | we have that, for any a G ©(2p), 

tt((7i72)"^a/i,2) = tl(( 7 i/i 72 / 2 )"^a) + 2fc- |Im(/i)| - |Im(/ 2 )|. 


We can thus derive from the graphical calculus for Wishart matrices (in complete analogy to the way formula (18 \ 
was obtained) that, for any d G N, 

\ p-i 3 2 


(40) 


E 


WAB'^Wj2,,d2 


Tr 


WAB*=(j)j 


/i:{l....,p}->-[fc] aee ( 2p ) 

/2:{p+l....,2p}->[fc] 


U/1./2) 
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where for each a £ &(2p) and /i ,p} —>■ [A], /2 : {p + 1,..., 2p} [fc], 

n{a,fij2) = 211(a) +jl(( 7 i 72 )"^a) + jl((7i/i 72 / 2 )"^a) +2fc- |Im(/i)| - |Im(/ 2 )|. 

Yet, by Lemma |A. 1| and equation (301 in Lemma | A.5 [ we get: First, for any a £ &{2p), 

(41) D(a) + D(( 7 i 72 )“^a) = 4p - (|a| + |( 7 i 72 )“^Q;|) ^ 4p - I 7172 I = 2p + 11 ( 7172 ) = 2p + 2, 

with equality if and only if a = aia 2 where a-i £ NC({1,... ,p}), ao £ JVC({p+ l,...,2n|). And second, for any 
a £ 6 ( 2 p) and any /i : { 1 ,... ,p} -)► [k], f 2 : {p + 1,... ,2p} ^ [k], 

(42) tl(a) + |l((7i/i72/2)"^a) < 2p + 11 ( 71 / 172 /J = 2p+ |Im(/i)| + |Im(/ 2 )|, 

with equality if and only if a = aia 2 where ai £ NC{{1,... ,p}) and /i o ai = /i, a 2 £ NC{{p + 1,..., 2p}) 

and /2 o a 2 = / 2 - So putting equations (411 and (421 together, we get in the end that for any a £ &(2p) and 

/i : { 1 , ■ ■ ■ ,P} [k], /2 : {p + 1 , ■ • ■, 2 p} ^ [k], 

n(a, fi,f 2 ) < 4p + 2fc + 2, 

with equality if and only if a = aia 2 where ai £ NC{{1,... ,p}) and /i o ai = /i, 02 S NC{{p + 1,..., 2p}) and 
/2 O a 2 = / 2 - 


We thus get that, asymptotically, the dominant term in formula (401 factorizes as 

2 


E Tr 


wab'= (j) 


d—>-+oo 


^4p-\-2k-\-2 


E 


E 


Jt(aia 2 ) 


QiGA''C({1,...,p}) /i:{l,...,p}-l[fc], /ioai=/i 

a2eNC({p+l,...,2p}) /2:{p+l,...,2p}->'[fe], /20Q2=/2 


/ 


d—>-+oo 


^2p+fc+l E E 


r.tt(“) 


afciVOlpj/:[PJ-1[/CJ 

\ /oo, = / / 

where the last equality is simply because ll(aia 2 ) = |l(ai) + tl(a 2 )- And hence, 

PIN 2 


EiVab^Wj 2 rf 2 I T'' 


E^AB^(j)j 


d—>-+oo 


^2p+fc+l E (cfc)**^“^ 

aeNC{p) 




which is exactly what we needed to conclude the proof. 


□ 


Let us illustrate the result stated in Theorem |H.1| in the simplest case of 2-extendibility and uniformly distributed 
mixed states. In Figure]^ the spectral distribution of Wd = (Wabi ® Wbs + Wab2 ® Wbi) /d'^, for Wab ~ yVd^,d'^, 
and a Marcenko-Pastur distribution of parameter 2 are plotted together. The empirical eigenvalue histogram is 
done in dimension d = 12 , from 100 repetitions. 

H.2. “Modified” GUE ensemble. Fix k £ N. Then, for each d S N, let G ~ GUE{d^) and define the random 
Hermitian matrix Gd on (0“^)®^+^ by 


(43) 






In complete analogy to what was explained in the case of Wishart matrices. Proposition |2.3| establishes that when 
d —> +CX), the eigenvalue distribution of Gd converges in moments towards a centered semicircular distribution of 
parameter k. But here again, there is in fact convergence in probability of towards psc{k)j which is made 
precise in Theorem H.3 below. 

Theorem H.3. For any I £ R and any e > 0, 

lim Pc-GUEid^) UNgAI) - PSG(k){I)\ > e) =0, 

Cl —^~1“00 


where the matrix Gd is as defined in equation 

As already explained in the Wishart case, this follows directly from the moment’s estimate in Proposition |2.3[ 
together with the variance’s estimate, for all p £ N, 

2p1 


(44) 


'^^GAB~GUE{d^) Tr 


|^E^AB'°(i) 


= o(d^P+^+A. 

d — 



















fc-EXTENDIBILITY OF HIGH-DIMENSIONAL BIPARTITE QUANTUM STATES 


33 


Figure 8. Spectral distribution of (T4 ^abi ® Wbs + ^abs ® WbJ for Was, ^ y^(p,cp vs 
Marcenko-Pastur distribution of parameter 2. 



The proof follows the exact same lines as the one of Proposition |H.2| and is not repeated here. 

Let us illustrate the result stated in Theorem in the simplest case of 2-extendibility. In Figure the 
spectral distribution of Gd = (GaBi ® Wbs + Gabs ® IdeJ /(P, for Gab GUE{P), and a centered semicircular 
distribution of parameter 2 are plotted together. The empirical eigenvalue histogram is done in dimension d = 10, 
from 100 repetitions. 


Figure 9. Spectral distribution of (Gabi ® IdBa + Gabs ® Wbi) for Gab ~ GUE(d^) vs 
Centered semicircular distribution of parameter 2. 
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Remark H.4. Wt may in fact say even more on the convergence of the random matrix sequences (W£;)(jgN O'^d 
(Gd)dgN defined by equations (391 and (43) respectively. Namely, 


^Wd , tJ-MP(ck) and Nod . —> h-sc(k)- 

d —^-j-OO d —^-|-oo 


To establish this almost sure convergence result, the only thing that has to be verified is that, for any p G N, the 
series of variances 


+ 00 


(45) 


^Var 




1 


dk+1 


Tr(Vk|) 


+ 00 


and Var 




1 


dk+1 


Tr 


(g?) 


are summable. Indeed, almost sure convergence will then automatically follow from a standard application of the 
Chebyshev inequality and the Borel-Cantelli lemma. And condition (45) actually holds, as a consequence of the fact 
that, for any p S N, 


Var 


fjk+l 


Tr«) 


= O {d and Var 


1 


dk+1 


Tr 


(g?) 


= o(d-'‘). 


H.3. Asymptotic freeness of certain Gaussian matrices. Let us fix a few definitions and notation. Given 
n G N and a classical probability space, we define the free probability space \M.n{L°°[^,P\),Tn\, where 

A4„(L°°[n, P]) is the set oinxn matrices with entries in L°°[Tl,P] and ipn{') = ETr(-)/n is the normalized trace 
function on A4n{L°°[^, P])- The two particular examples we shall focus on in the sequel are the ones we have 
already been extensively dealing with, namely GUE and Wishart matrices. 

Lemma H.5. Given two finite-dimensional Hilbert spaces A = B = and G a random GUE matrix on 

A G B, we define the following random matrices on A G Bi ® B 2 : 


Gi = 


1 


GaBi G Wbj and G 2 = 


1 


y/d+A 


Gab2 G Wbi- 


Then, for any p G N and any function f : [2p\ -+ [2], 

(46) +00 (^/(i) • ■ ■ ^/(2p)) = (^ e NC^'^\ 2 p) : / o A = /| . 

Proof. We know from the proof of Proposition 2.3 (and using the same notation as those employed there) that 

^ \ ,)t(7-G)-p-l ,#(77A)-p-|Im(/)| 

TdAdl [Uf(l) ■ ■ ■ G/(2p) J - 2^ «A “b 

AGqjC)(2p) 

Now, as explained there as well, for any A G fp*'^^(2p), on the one hand j)(7“^A) ^ p + 1 with equality iff A G 
NC^'^\2p), and on the other hand {((y^^A) ^ p + | Ini(/)| with equality iff A G NCdU{2p) and / o A = /. The 
asymptotic estimate (46) therefore immediately follows. □ 

Theorem H.6. Let Gab a random GUE matrix on A G B. Then, the random matrices Gabi G IdB 2 
Gab 2 G Wbi on a (g) Bi G B 2 are asymptotically free. 

Proof. Theorem H. 6 | is a direct consequence of Lemma H.5 (see e.g. [30], proof of Proposition 22.22, for an entirely 
analogous argument). Indeed, as d^, dn + 00 , the two empirical spectral distributions and pg^ both converge 
to the semicircular distribution with mean 0 and variance 1. And equation (46) is exactly the rule for computing 
mixed moments in two free such semicircular distributions (see e.g. EHI, Lecture 12). □ 

Lemma H.7. Given two finite-dimensional Hilbert spaces A = B = and W a random Wishart matrix 

on A g B with parameter cd^d^ G N, we define the following random matrices on A g Bi 


Bn; 


= 


1 


d\d^ 


ITabi G IdB2 and W 2 = 


1 


dAds 


1Tab2 G Wbi ■ 


Then, for any p G N and any function / : [p] —>■ [2], 


(47) 




J{<y) 




aGNC{p) 

/oa = / 
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Proof. We know from the proof of Proposition 6.2 (and using the same notation as those employed there) that 

aeG(p) 


Now, as explained there as well, for any a G @{p), on the one hand tl(a) + tl( 7 “^a) ^ p + 1 with equality iff 
a G NC{p), and on the other hand (((a) -|- 'i{'yj^a) < p + | Im(/)| with equality iff a G NC{p) and f o a = f. The 
asymptotic estimate (471 therefore immediately follows. □ 


Theorem H.8. Let Wab be a random Wishart matrix on A 0 B with parameter cdA^B G N. Then, the random 
matrices Wabi ® IdB2 Wab2 ® Wbi on A 0 Bi 0 B2 are asymptotically free. 


Proof. Theorem H.8| is a direct consequence of Lemma H.7 (see e.g. [3D], proof of Proposition 22.22, for an entirely 
analogous argument). Indeed, as dx, ^b +00, the two empirical spectral distributions and both converge 
to the Marcenko-Pastur distribution with parameter c. And equation (471 is exactly the rule for computing mixed 
moments in two free such Marcenko-Pastur distributions (see e.g. [30], Lectures 12 and 13). □ 
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